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droplet growth rate due to water vapor equilibrium, 
ym sec" 1 

droplet growth "kernel 11 for H2SO4-H2O condensation, 
ym cm 3 sec -1 molecule -1 

altitude grid spacing, 2xl0 5 cm 

discrete analog of particle size 

discrete analog of altitude 

coagulation kernel for particles of radii r and r f , 
cm 3 particle -1 sec -1 

Knuds en number 

coagulation kernel for particles of sizes i and j at 
altitude level k, cm 3 particle -1 sec -1 

discrete analog of time 

sulfuric acid vapor loss rate due to condensation on par- 
ticles, sec* 1 

molecular mass of species v, g 
molecular weight of species v, amu 

aerosol droplet size distribution, particles cm -3 ym -1 

boundary concentration for species v, molecules cm* 3 

concentration of aerosol droplets of size i, particles 
cm -3 


1 In the text some units have been modified for consistency in certain 
expressions . 
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total concentration of air molecules, molecules cm -3 

concentration of the gaseous species v , molecules cm” 3 
(v = a for H2SO4, v = H for H 2 0 ) 

equilibrium vapor concentration (vapor pressure) of 
species v over an aerosol droplet, molecules cm” 3 

condensation nuclei size distribution, particles cm -3 pm' 


concentration of condensation nuclei of size i, 
particles cm” 3 

Avogadro ? s number (6.02xl0 23 molecules mole -1 ) 

condensation nuclei production rate from evaporating 
droplets, particles cm -3 ym” 1 sec” 1 

sulfuric acid vapor production rate, molecules cm -3 sec” 1 

particle radius, ym 

heterogeneous heteromolecular nucleation rate, sec” 1 
aerosol droplet rainout rate, sec -1 
time, sec 

air temperature, K 

total core volume of the droplets n(r,z,t), ym 3 cm” 3 ym” 
total core volume of the droplets n^ ? ym 3 cm~ 3 

boundary recombination velocity, cm sec -1 

particle gravitational sedimentation velocity, cm sec -1 

particle volume, ym 3 

total particle volume per unit volume of air at height k 
and time £ , ym 3 cm” 3 

weight fraction of sulfuric acid in aqueous solution 
altitude, cm 

sticking coefficient for H2SO4 molecules on aerosol 
droplets 
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total core volume second moment of the droplets 
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size distribution of droplets according to their 

volume, v, and their core volume w, particles cm” 3 yn ” 6 

double counting factor for coagulation terms (0=1, i ^ j; 
1/2, i = j) 

aerosol stratification factor 

aqueous sulfuric acid solution density, g cm” 3 
aerosol droplet core density, g cm” 3 
average aerosol particle density, g cm” 3 
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numerical time step, sec 
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growth time of droplets from size i to size i + 1, sec 

water vapor growth time for droplets from size i to size 
i + 1, sec 




vertical particle diffusion flux, particles cm 2 ym 1 sec 1 


vertical particle sedimentation flux, particles 

—o —1 -1 

cm * ym sec - 1 



total vertical particle flux due to sedimentation and 
diffusion, particles cm” 2 ym” 1 sec 1 
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vertical droplet core volume flux, ym 3 cm 2 ym 1 sec 1 
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SUMMARY 


A time- dependent one-dimensional model of the stratospheric sulfate 
aerosol layer is presented. In constructing the model, a wide range of basic 
physical and chemical processes are incorporated in order to avoid predeter- 
mining or biasing the model predictions. The simulation, which extends from 
the surface to an altitude of 58 km, includes the troposphere as a source of 
gases and condensation nuclei and as a sink for aerosol droplets. The size 
distribution of aerosol particles is resolved into 25 categories with particle 
radii increasing geometrically from 0.01 to 2.56 pm such that particle volume 
doubles between categories. In the model, sulfur gases reaching the strato- 
sphere are oxidized by a series of photochemical reactions into sulfuric acid 
vapor. At certain heights this results in a supersaturated H2SO4-H2O gas mix- 
ture with the consequent deposition of aqueous sulfuric acid solution on the 
surfaces of condensation nuclei. The newly formed solution droplets grow by 
heteromolecular heterogeneous condensation of acid and water vapors; the 
droplets also undergo Brownian coagulation, settle under the influence of 
gravity, and diffuse in the vertical direction. Below the tropopause, parti- 
cles are washed from the air by rainfall. A set of continuity equations is 
derived which describes the temporal and spatial variations of aerosol drop- 
lets and condensation nuclei concentrations in air as well as the sizes of 
cores in droplets; numerical techniques to solve these equations accurately 
and efficiently are also developed. Calculations illustrate the behavior and 
precision of the computer algorithms. The model correctly predicts many of 
the observed properties of the sulfate layer, including its height, density, 
composition, and particle size spectra. 


INTRODUCTION 


The stratospheric sulfate layer is a highly dispersed aerosol of sulfuric 
acid-water solution droplets which is usually found between altitudes of 15 
and 30 km. The layer was discovered by Junge and his coworkers nearly two 
decades ago (ref. 1) and since then has been a subject of considerable inter- 
est to climatologists and atmospheric scientists. This interest has been 
stimulated by recent suggestions that the aerosol may exert an important 

*Now at Systems and Applied Sciences Corporation, Hampton, Virginia 23666 . 
tNow at Georgia Institute of Technology, Atlanta, Georgia 30332 . 
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influence on the climate of Earth (refs. 2-4) as well as being responsible for 
brilliant sunsets and some rather unusual optical phenomena (ref. 5). Addi- 
tionally, there is concern that anthropogenic perturbations of the layer due 
to increased sulfur loading of the stratosphere or the injection of particu- 
late matter at high altitudes could possibly affect the climate or catalyze 
heterogeneous chemical interactions that affect global ozone abundances 
(ref. 6). As a consequence, there has been a fairly concerted effort to study 
the layer, both experimentally and theoretically, and to develop analytical 
tools to evaluate the effect of perturbations on the layer. In this paper we 
describe a time-dependent one-dimensional model that accurately simulates many 
of the observed properties of the stratospheric aerosol layer and may be 
useful as a predictive device. 

Because of the large number of observational studies that have been car- 
ried out since its discovery, the characteristics of the sulfate layer are 
relatively well-defined. Many investigators have used balloon-borne diffusion 
chambers (refs. 1 and 7) and optical particle counters (refs. 8-11) to analyze 
the aerosol particles. Other experiments have utilized high-flying aircraft 
equipped with wire impactors or filters to collect particle samples 
(refs. 12-15). Sunlight scattered by the layer has been detected from satel- 
lites (ref. 16) and the layer has been extensively probed with laser light 
using lidar techniques (refs. 17 and 18). The results of these experiments 
have led to a quite comprehensive and consistent picture of the aerosol layer. 

In the sulfate layer, the height of the maximum in the number mixing 
ratio of large aerosol particles (radii greater than about 0.15 ym) usually 
lies between 20 and 25 km and appears to be closely related to the altitude 
of the tropopause (ref. 8). The aerosol mass mixing ratio is also peaked 
between 20 and 25 km in the stratosphere. The sulfate layer exhibits sea- 
sonal, latitudinal, and short-time scale variability and can be greatly 
enhanced following major volcanic activity (ref. 19). Aerosol particle size 
distributions fall off rapidly for particle radii greater than 0.1 ym, and 
the ratio of the concentrations of particles with radii greater than 0.15 ym 
to those greater than 0.25 ym is about three to five in the layer. Unlike 
the mixing ratio of large particles, the number mixing ratio of total parti- 
cles falls off rapidly above the tropopause. 

A model of the stratospheric sulfate aerosol should reproduce the general 
characteristics of the aerosol layer indicated above. In the past, aerosol 
models have, for the most part, been developed to investigate specific aerosol 
properties. Junge et al. (ref. 1) first used a one-dimensional diffusion 
model with particle sedimentation and coagulation to explain the observed 
rapid decrease of the total aerosol particle mixing ratio above the 
tropopause. Friedlander (ref. 20) formulated a detailed model of particle 
coagulation and applied it to the prediction of stratospheric aerosol size 
distributions. Later, Martell (ref. 21) utilized a model similar to Junge’ s 
to study the transport and coagulation of Aitken (small) particles in the 
lower atmosphere. 

Over the past decade, a large number of aerosol models dealing princi- 
pally with aerosol transport has been published. For example, Hunten 
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(ref. 22) used a one-dimensional diffusion-sedimentation model for a monodis- 
persed aerosol to estimate the residence times of particles in the strato- 
sphere. Hoffmann et al. (ref. 23) employed a similar model, but for several 
different particle sizes, to calculate the stratospheric accumulation of alu- 
minum oxide particles emitted by Space Shuttle booster engines. Cadle et al. 
(refs. 24 and 25) investigated the global dispersion of aerosol particles fol- 
lowing volcanic activity using a two-dimensional model with particle sedimen- 
tation and parameterized meridional transport. 

Some recent aerosol models have treated aerosol physics in more detail. 
Kritz (ref. 26) presented a novel and detailed treatment of aerosol growth and 
coagulation and investigated the development of stratospheric particle size 
distributions by these processes. Burgmeier and Blifford (ref. 27) con- 
structed a "box" model for the stratospheric aerosol, which includes particle 
coagulation and growth (by absorption of sulfur dioxide molecules) . They used 
their model to study aerosol "aging" — that is, the short-time evolution of 
different initial particle size distributions. Rosen and Hofmann (refs. 28 
and 29) extended their earlier work (ref. 23) to account for particle coagula- 
tion and growth, and they reestimated the effect of Space Shuttle operations 
on the stratospheric aerosol layer. 

Photochemical theories involving atmospheric sulfur constituents have 
developed in parallel with the physical aerosol models. Junge et al. (ref. 1) 
first suggested that the large aerosol particles observed in the stratosphere 
might be formed by the condensation and reaction of vaporous sulfur compounds. 
Subsequently, Scott et al. (ref. 30) and Friend et al. (ref. 31), among 
others, introduced chemical reaction mechanisms by which sulfur dioxide mole- 
cules absorbed by aerosol droplets could be oxidized into sulfate. Other 
workers have proposed gas phase reaction sequences that can convert sulfur 
oxides into sulfuric acid vapor (e.g., refs. 15, 32, and 33) which readily 
condenses on aerosol particles (e.g., ref. 34). Using one-dimensional trans- 
port models, Harrison and Larson (ref. 35), Harker (ref. 36), and Crutzen 
(ref. 33) have all studied atmospheric gas phase sulfur chemistry and have 
found that sulfate can be produced in the lower stratosphere at a rate fast 
enough to maintain the observed aerosol layer. 

In the remainder of this report, we construct a new time-dependent one- 
dimensional model of the stratospheric sulfate layer. In our model, we have 
attempted to combine all of the important features of aerosol physics and 
sulfur photochemistry found in earlier models. In addition, we have adopted 
the most recent theories for the nucleation and growth of aerosol droplets as 
well as some newer measurements of the abundances and reactivities of sulfur 
compounds in air. We have also introduced a treatment of aerosol cores that 
describes their accumulation and evolution in stratospheric solution droplets; 
the analysis of core size and composition can be a valuable theoretical and 
diagnostic tool for understanding the aerosol layer. Our model is a compre- 
hensive, interactive gas-particle simulation of the stratospheric aerosol. 

In the second section of this paper. Physical Chemistry of Atmospheric 
Sulfur, we consider the sources and sinks of atmospheric sulfur compounds, 
their photochemical interactions in air, and the relation of gaseous sulfur 
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chemistry to the formation of sulfuric acid solution droplets. The physical 
aerosol processes that are incorporated into our model are discussed in the 
third section, Physical Aerosol Processes, and some related parameter values 
are tabulated in appendix A. The continuity equations for aerosol particles 
and condensation nuclei and for the cores within the aerosol droplets are 
derived in the fourth section. The Aerosol Continuity Equation, and in appen- 
dix B. Our treatment of coagulation and of particle growth by condensation is 
also outlined in the third and fourth sections and in appendix B. In the 
fifth section, A Numerical Representation of the Aerosol Continuity Equations, 
and appendix C, we establish the finite difference equations in one (vertical) 
direction and solve them to obtain time-dependent aerosol particle distribu- 
tions. In the sixth section. Model Validations and Simulations, we present 
the results of model validation tests and briefly compare our model predic- 
tions of aerosol properties with stratospheric observations. Tables 2-4 sum- 
marize most of the important sources of uncertainty (chemical, physical, and 
numerical) for our model calculations. Finally, conclusions about the poten- 
tial applications of our aerosol model are given in the seventh section. 

In a companion publication (ref. 37) hereafter referred to as NASA 
TP-1363, we make a more detailed comparison of our model predictions with 
aerosol observational data. In NASA TP-1363, we also present a comprehensive 
sensitivity analysis showing areas where our aerosol model is most uncertain 
and suggest new observations that might help to resolve important questions 
about the origin of the stratospheric aerosol layer. 


PHYSICAL CHEMISTRY OF ATMOSPHERIC SULFUR 


Gaseous sulfur compounds have an important role in the formation and 
maintenance of the stratospheric aerosol layer. One goal of our aerosol model 
is to investigate the interactions of sulfur gases with aerosol particles. 
Accordingly, we need to determine the abundances and distributions of volatile 
sulfur constituents in the environment. In this section we discuss the known 
sources, sinks, and concentrations of atmospheric sulfur-bearing molecules and 
outline the stratospheric photochemistry applicable to these species. The 
aeronomical model, which we use to compute the concentrations of sulfur com- 
pounds in air, is also described. Finally, we present and discuss typical 
vertical distributions of the sulfur constituents as predicted by our model. 


Sulfur Compounds in the Atmosphere 

Sulfur gases can be introduced into the atmosphere in several different 
forms by a host of sources. For example, industry releases copious amounts of 
sulfur dioxide (SO 2 ), and smaller amounts of hydrogen sulfide (H 2 S), carbon 
disulfide (CS 2 ), and carbonyl sulfide (OCS) into the environment 
(refs. 38-40). Sulfur oxides, including sulfuric acid vapor, can be generated 
in automotive engine exhaust systems. Hydrogen sulfide is a by-product of the 
decay of vegetable matter (ref. 41). Dimethyl sulfide, (CH 3 ) 2 S, may be pro- 
duced in the Earth* s oceans (refs. 42 and 43), and carbon disulfide has been 
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found dissolved there (ref. 44). Volcanic eruptions also inject sulfurous 
compounds into the air (e.g., ref. 45). These and a few other natural and 
anthropogenic sources probably account for most of the sulfur input to the 
atmosphere (ref. 41). 

Sulfur dioxide and hydrogen sulfide molecules can react in air with 0, 

0 2 , O 3 , and OH (ref. 46), leading to H 2 SO 4 and other stable sulfates (ref. 47). 
At low altitudes, water soluble gases like SO 2 , H 2 S, and H 2 SO 4 can be washed 
out of the atmosphere by rainfall (refs. 48 and 49); SO 2 and H 2 S absorbed in 
hydrometeors can react in solution with dissolved O 2 and O 3 as well as with 
other substances to form soluble and insoluble sulfated compounds (refs. 38, 
47, and 50). Near the ground, sulfur gases can be removed on plant and min- 
eral surfaces or be dissolved in the seas (refs. 41 and 47). The sulfur 
compounds containing H-C bonds are probably attached by hydroxyl radicals in 
the troposphere, eventually forming soluble sulfur oxides. 

Sulfur dioxide and hydrogen sulfide gases have been detected at low alti- 
tudes, but they exhibit a wide range of concentrations with sea-level values 
usually less than a few parts per billion by volume (ppbv) and often very 
much less (refs. 41, 47, 51, and 52); and the effects of meteorology, urban 
pollution, and gas analysis techniques on these measured concentrations are 
unknown, however. 

Based on sulfur mass balance calculations, the tropospheric residence 
time for SO 2 has been estimated to be only a few days (refs. 41 and 53). 
Hydrogen sulfide, because of its extreme reactivity, is estimated to have an 
atmospheric lifetime much less than a day (refs. 47 and 54). Likewise, sul- 
furic acid vapor and dimethyl sulfide are also expected to have very short 
low altitude residence times. Nevertheless, Jaeschke et al. (ref. 55) have 
observed about 0.05 ppbv of SO 2 in the upper troposphere and lower strato- 
sphere. The altitude-dependence of the tropospheric SO 2 mixing ratio observed 
in reference 55 appears to be consistent with the estimates by Friend 
(ref. 41) and Levy (ref. 53) of the tropospheric lifetime of SO 2 . 

Carbonyl sulfide has been measured near the ground over a wide geographi- 
cal range. Hanst et al. (ref. 56) found about 0.2-0. 3 ppbv of OCS in air; 

A. Bandy (private communication, 1977) detected about 0.47 ± 0.05 ppbv; and 
Sandalls and Penkett (ref. 57) observed about 0.51 ± 0.05 ppbv. Disregarding 
the Hanst measurement, which has a large experimental uncertainty associated 
with it, OCS appears to have a very uniform concentration of about 0.5 ppbv 
throughout the lower troposphere (according to the authors of ref. 57, a large 
fraction of their reported OCS variability may actually have resulted from 
their gas analysis procedure). Unlike most other sulfur gases, OCS is 
expected to have a long tropospheric lifetime, possibly a year or more. 

Carbon disulfide is a less inert chemical than OCS, but it may still have 
some aeronomical significance because of its possible oxidation to OCS 
(ref. 58). Sandalls and Penkett (ref. 57) measured an average CS 2 concentra- 
tion of about 0.10 ppbv near the surface, with a very large variability 
(±50%). 
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By contrast with observations in the lower troposphere, indications are 
that, except for OCS and SO 2 , sulfur molecule concentrations in the upper 
troposphere are very low (refs. 55 and 59). 

Because of the apparent instability of most of the dominant sulfur com- 
pounds in air, and the corresponding failure to detect large amounts of them 
at high altitudes, it seems that only a small fraction of the total quantity 
of these substances which are emitted at the ground actually reach the strato- 
sphere and thus influence the aerosol layer. Accordingly, in the present 
model, we consider only S0 2 and OCS as gaseous sources of sulfur for aerosol 
particles. 

To simulate the small percentage of the surface emission of S0 2 that is 
transported to the lower stratosphere, we have adjusted its ground mixing 
ratio to about 0.6 ppbv and allowed the S0 2 molecules to diffuse upward 
through the tropopause. After depletion by rainout and chemical oxidation, 
the S 0 2 abundance near the tropopause is about 0.05 ppbv, in agreement with 
the high-altitude measurement of reference 55. There are several mechanisms 
by which sulfur dioxide may be directly injected into the stratosphere, 
including deposition by aircraft engines and volcanic eruptions (ref. 45). 
However, for the purpose of establishing the ambient properties of the aerosol 
layer, we have not included sulfur gas injection in our present calculations. 

As we have already mentioned, OCS originates in biological, volcanic, and 
industrial processes (refs. 40 and 60) and is transported upward into the 
stratosphere, much like the halocarbons CFCI 3 and CF 2 C1 2 (ref. 61), where it 
can be photocheraically transformed into S0 2 in situ (ref. 33). To determine 
the quantity of sulfur dioxide that is generated by OCS photolysis, we have 
adopted a tropospheric OCS concentration of 0.47 ppbv and the OCS transport 
and photolysis rates to be discussed shortly. 

The total weight of sulfate leaving the stratosphere in aerosol has been 
estimated by Lazrus and Gandrud (ref. 14) to be equivalent to 1.2x10^ 
metric tons of S0 2 per year and by Hofmann et al. (ref. 8) to be about 
2.5 times greater (1x10 s metric tons of S0 2 annually is equal to about 
6.5xl0 6 S0 2 molecules per square centimeter column per second, globally aver- 
aged). The mass of S0 2 that diffuses into the stratosphere through the tropo- 
pause in our model is almost 2xl0 5 metric tons yearly for our assumed S0 2 
concentrations and vertical diffusion coefficients (which will be discussed 
later in this section). However, most of this S0 2 is converted into particle 
sulfate just above the tropopause level and is removed rapidly from the 
stratosphere as particulate matter. At a distance a few kilometers above the 
tropopause, the upward diffusive flux of tropospheric sulfur dioxide is 
greatly reduced, and the downward flux of aerosol sulfate is due primarily to 
OCS decomposition at higher altitudes (equivalent to about 9x10^ metric tons 
of S0 2 per year). 


The Photochemistry of Sulfur Compounds 

Sulfur dioxide, carbonyl sulfide, and other sulfur gases that reach the 
stratosphere can be photochemically oxidized into sulfuric acid vapor. 
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Table 1 presents a set of photochemical reactions by which OCS molecules are 
transformed into H 2 SO 4 molecules; these are the sulfur reactions used in our 
model. 

Since OCS is a very stable compound, its aeronomy is quite limited. In 
fact, besides its source at the ground, only photolysis and transport are 
thought to influence its distribution in the upper atmosphere (ref. 33). The 
OCS photoabsorption coefficients have been measured (refs. 70 and 71). The 
observed absorption spectrum, illustrated in figure 1 , when combined with the 
solar fluxes of Ackerman (ref. 76), leads to the daily averaged OCS photodis- 
sociation rates given in table 1 . 

According to reference 77, the OCS ultraviolet photolytic process is 

OCS + hv + S + CO (1) 

where the sulfur atom can be electronically excited. We do not distinguish 
between the different states of the product sulfur atom here since its degree 
of excitation does not affect our results. 

Carbonyl sulfide reacts very slowly with oxygen atoms in the atmosphere 
(see reaction RIO in table 1). In our model, less than 10% of the OCS reach- 
ing the stratosphere is decomposed by this reaction. The reaction of OCS with 
excited oxygen atoms, 0( 1 D), is even less important; although the rate coeffi- 
cient is probably large (~10 -10 cm 3 sec -1 ), 0( 1 D) concentrations are extremely 
small in the stratosphere (<100 atoms cm -3 ). The reaction of OCS with OH is, 
likewise, negligible; in a recent experiment (ref. 78), an upper limit of 
7xlO “ 15 cm 3 sec -1 was placed on the room temperature rate constant. We have 
not considered other possible (tropospheric) loss processes for OCS in this 
study because our aim here is simply to determine the stratospheric distribu- 
tion of OCS when its tropospheric mixing fraction is known. 

Carbon disulfide will be decomposed in the atmosphere by processes analo- 
gous to, but much faster than, those for OCS (e.g., ref. 40). Tropospheric 
CS 2 concentrations are lower and more uncertain than those of OCS. Moreover, 
the fraction of CS 2 reaching the stratosphere is probably much less than that 
of OCS. Accordingly, in this report we have omitted CS 2 as an active strato- 
spheric sulfur constituent, recognizing that its role in forming the sulfate 
aerosol layer may need to be evaluated later. 

The sulfur atoms generated during OCS photodecomposition react rapidly 
with molecular oxygen to form SO. Sulfur monoxide is, in turn, rapidly oxi- 
dized into SO 2 . There are three reactions (R5, R 6 , and R7) that further oxi- 
dize SO 2 * and each has been proposed at one time or another as an important 
mechanism leading to stratospheric H 2 SO 4 production. Except for reaction (R7) 
these processes are all quite slow, having time constants ranging from months 
to years; the rate of S0 2 oxidation by reaction (R7) is dominant by at least 
one order of magnitude throughout the lower stratosphere. 

The sulfur radical HSO 3 formed by reaction (R7) can react with 0 2 to pro- 
duce HSO 5 , which can react in turn with nitrogen oxides to yield compounds 
such as HOSO 2 ONO and HOSO 2 ONO 2 (ref. 79). Because the photochemistry of these 





latter species is not well known, we have not been able to include them in our 
model. However, it is interesting to note that Farlow et al. (ref. 80 ) have 
recently detected small amounts of similar nitrogen-bearing sulfate compounds 
in stratospheric aerosol samples. 

As a radical in air, HS 0 3 is likely to be quickly transformed into a more 
stable H2SO4 molecular configuration. Therefore, it is reasonable to assume 
that the sequence of reactions leading from HSO3 to H2SO4 does not limit the 
overall rate of conversion of SO2 into sulfuric acid vapor. We have proposed 
reaction (R8) — which may actually occur in the atmosphere, although at a 
slower rate than the direct oxidation of HSO3 by O2 — followed by reaction 
(R 9 ) as a substitute for a more complex HS03-to-H2SC>4 reaction sequence (also 
see ref. 81 ). In fact, detailed reaction schemes have been suggested; some 
result in H2SO4 formation and some do not, but none of them is supported by 
rate constant data (refs. 32 and 82 ). Graedel (ref. 83 ) has considered sev- 
eral HSO radicals in an analysis of the role of sulfur oxides in photochemi- 
cal smog and has roughly estimated rate constants for a hypothetical series of 
HSO reactions. However, regarding the importance of these species to aerosol 
formation, Graedel’ s work is not conclusive and, in addition, his proposed 
reaction scheme cannot be realistically applied under stratospheric condi- 
tions. For these reasons we have chosen to defer consideration of any addi- 
tional sulfur radical reactions until more chemical rate information is 
available. 

In our model, all of the reactions we have discussed so far lead to the 
creation of sulfur trioxide, SO3 , which in turn rapidly combines with H 2 0 to 
form sulfuric acid vapor (reaction R 9 ) . Considering that stratospheric water 
vapor concentrations are about 3-6 ppmv (ref. 84 ), it is doubtful that SO3 
has any other significant aeronomical interactions. For example, its photol- 
ysis rate would have to exceed 0.1 sec -1 to compete with reaction (R 9 ) any- 
where within the aerosol layer. There is still some doubt about the mechanism 
of the reaction between SO3 and H2O; it may be a heterogeneously catalyzed 
process. Accordingly, the role of S 0 3 as a precursor of sulfuric acid vapor 
might be exaggerated in our reaction scheme. Until an alternative gas-phase 
sulfur oxidation mechanism is verified experimentally, however, the reactions 
in table 1 are sufficient for our purposes. 

Stratospheric sulfuric acid molecules generated by SO3 hydration can dif- 
fuse to high altitudes and be photodissociated , presumably forming SO2 in one 
or more photochemical steps; they can diffuse to the troposphere and be fil- 
tered from the air by rainfall, or they can interact in situ with aerosol 
particles. H2SO4 absorption coefficients in the ultraviolet spectrum are 
unavailable. We have arbitrarily adopted HC 1 absorption cross sections 
(ref. 75 ) for H2SO4 even though H2SO4 binds hydrogen atoms more tightly than 
HC 1 . Obviously, a measurement of the H2SO4 ultraviolet absorption coefficient 
is needed. But even so, the photolytic lifetime of sulfuric acid vapor below 
50 km is probably long enough (several months or more) that chemical, aero- 
sol, and transport processes can control its concentration in the strato- 
sphere. Photodissociation will be more important at higher altitudes. In 
this case, when H2SO4 is being fragmented by energetic photons, it is prob- 
ably a reasonable approximation (ignoring other possible dissociation prod- 
ucts) to assume that the initial step in the photodecomposition releases S 0 2 . 
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In the mesosphere, where the radiation environment is more intense, the SO2 
can be further reduced to SO or S. 

The gas phase reactions of H2SO4 molecules with the atoms and radicals 
normally found in the atmosphere have not been studied experimentally. 

Although it is possible that H2SO4 could react with 0 , Cl, and OH in the 
middle and upper stratospheres, no rate constant data are available for these 
interactions. The H2SO4 photolysis rate is quite small throughout the strato- 
sphere (see table 1 for the values corresponding to our assumed H2SO4 disso- 
ciation cross sections). Accordingly, even a very slow reaction — most 
probably one with OH — could affect the H2SO4 abundance in the upper strato- 
sphere. On the other hand, such a reaction would not materially affect our 
predictions for the aerosol layer because at the higher altitudes where a 
slow radical reaction might be important, the aerosol layer is already quite 
evanescent. 

Molecular rainout and washout processes are difficult either to describe 
or model realistically. Tropospheric hydrometeor ensembles have a signifi- 
cant function as filters of sulfur gases migrating from surface sources to 
the stratosphere. In this work, we have, by necessity, only partially treated 
the problem of the tropospheric sulfate balance and the possible leakage of 
biospheric sulfur to the stratosphere. The usual aeronomic approximation for 
the precipitation removal of a gaseous species is to specify a decreasing 
rainout rate with increasing altitude such that the resulting average atmo- 
spheric residence time corresponds roughly to one deduced using global mass 
budget arguments. For SO2 we have adopted a lifetime against washout (plus 
rainout) of about 3 days at the surface; the rate of SO4 washout (or the 
inverse of the SO2 washout lifetime) is assumed to decrease linearly with 
increasing height above sea level to 0 at 13 km (see table 1 ). Presently, 
there are insufficient data to evaluate the tropospheric residence time of 
H2SO4 against heterogeneous removal; we have therefore assigned it a washout 
lifetime of 0.5 day at ground level and a washout profile similar to that 
for SO2 - The radical HS 0 3 is assumed to be washed from the air at the same 
rate as H2SO4. (Later, we will consider tropospheric precipitation as a sink 
for downward diffusing aerosol droplets as well.) 

The interaction of sulfuric acid vapor with stratospheric aerosol drop- 
lets is discussed in detail in the next section. H2SO4 condensation on and 
evaporation from aerosol particles are taken into account in our photochemical 
model. However, we have not included the heterogeneous interaction of other 
gaseous species with aerosol particles. Friend et al. (ref. 31 ) have pointed 
out that SO2 could react directly with sulfuric acid-water droplets to produce 
additional dissolved sulfates; the process would require the presence in solu- 
tion of a buffering agent such as ammonium ions. Recently, Eggleton and Pen- 
kett (ref. 85 ) observed that SO2 can be readily oxidized in aqueous solution 
by hydrogen peroxide and to a lesser degree by ozone; the rate of oxidation 
is sensitive to the solution composition, however. 

The inclusion of heterogeneous SO2 oxidation mechanisms in a model would 
obviously increase the rate at which sulfur gases are converted into aerosol 
particles. In NASA TP- 1363 , however, we have varied this conversion rate to 
show that it does not greatly affect the aerosol layer. 
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HS0 3 and other sulfur radicals might also condense or react on aerosol 
surfaces before they can be oxidized into H 2 SO 4 by gas phase reactions or 
might nucleate to form new aerosols. Unfortunately, nothing is presently 
known about the interaction of these species with stratospheric particles. In 
a test, we allowed HSO 3 to condense on aerosol particles at the same rate as 
H 2 SO 4 , and we assumed that HSO 3 reacts instantaneously in solution to form 
H 2 SO 4 . As expected, the concentrations of both HSO 3 and H 2 SO 4 (whose precur- 
sor is HS0 3 ) are reduced. However, the observable characteristics of the 
aerosol layer are essentially unaltered (neglecting possible effects on aero- 
sol nucleation mechanisms) . 

According to the previous discussion, it is reasonable to omit a detailed 
treatment of heterogeneous interactions from our present aerosol model while 
more suitable theories describing these interactions are being developed. All 
of the processes involving gaseous sulfur compounds that are included in our 
model are summarized in figure 2 . 


The Aeronomical Model 

We have used a one-dimensional model to calculate the atmospheric concen- 
trations of the sulfur constituents S, SO, SO 2 , SO 3 , HSO 3 , H 2 SO 4 , and OCS. 

The model includes the photoreaction scheme in table 1, augmented by the 
interaction of sulfuric acid vapor with aerosol particles. We have assumed 
that the vertical transport of gaseous species occurs by- eddy diffusion pro- 
cesses. Our technique for solving the related time-dependent one-dimensional 
species continuity equations is identical to a method we developed for a much 
larger aeronomical system (e.g., refs. 86-90). Summarizing our technique, we 
solve the species continuity equations at 30 equally spaced grid points 
between 0 and 58 km with a numerical finite difference formulation that is 
fully implicit in the transport terms and partially implicit in the photo- 
chemical terms for greater stability. (Implicit, in the numerical sense, 
means evaluation at the end of a computational time step. Implicit varia- 
bles are therefore quantities to be calculated during that step. Explicit 
variables are evaluated at the beginning of the step and are known) . 

We define a family of sulfur constituents, denoted as S0 X , which consists 
of all the sulfur species listed above except OCS. The photochemical reac- 
tions that affect SO include only those processes which give a net production 
or loss of sulfur in the S0 X family. Using the net SO chemistry, we compute 
the concentration of SO independently of the individual species concentra- 
tions. The S0 X abundance is then compared to the summation of the family 
species concentrations to estimate the cumulative error in the latter solu- 
tions and is used to correct those solutions if necessary (see ref. 87 for 
more details of this scheme). 

In our model, we also average the rate of each photochemical process 
diurnally to account for the effects of nighttime reactions on species abun- 
dances (ref. 90). In most respects, the aeronomical model used in this study 
is a simplified version of our more complete atmospheric model. 
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We have simulated vertical atmospheric transport in our model by the use 
of an "eddy" diffusion process. Our assumed eddy diffusion coefficient, shown 
in figure 3 and discussed in detail in reference 89, results in a very stable 
air layer near the tropopause which acts as a narrow escape valve for subsid- 
ing stratospheric gases and as an efficient barrier to rising tropospheric 
gases. Since the photochemical lifetime of S0 2 is not very long in the sul- 
fate layer, we expect that S0 2 injected or produced above the tropopause will 
not be extensively redistributed before it reacts and will almost always be 
chemically consumed near the altitude where it is added. 

We have selected the diffusion profile shown in figure 3 because refer- 
ence 91 has shown that it gives an excellent representation of the observed 
time-decay of excess stratospheric carbon-14 following atmospheric nuclear 
bomb tests. Our eddy diffusion coefficient also satisfies a second criterion: 
in an aeronomical model it leads to methane and nitrous oxide molecule con- 
centrations in the upper stratosphere that are close to measured values 
(e.g., refs. 89 and 92). Nonetheless, since there are large variabilities and 
uncertainties in atmospheric transport rates, the diffusion coefficients 
adopted here are certainly not unique. 

The actual mechanisms by which tropospheric gases and particles are 
transported into the stratosphere may be more closely related to injection 
than to diffusion. Such injection can occur — in an equivalent sense — during 
thundercloud penetration, tropopause folding, and vertical convection at the 
intertropical convergence zone. Horizontal advection would subsequently dis- 
tribute the injected material zonally and meridionally. Obviously, we cannot 
simulate the hydrodynamics of multidimensional transport with a one- 
dimensional model. However, we can roughly simulate the globally averaged 
mass transfer rates between the stratosphere and troposphere. More detailed 
studies of the meteorology of the stratospheric aerosol layer will have to 
wait the development of two- and three-dimensional dynamic models that can 
also accurately treat gaseous photochemistry and particle microphysics. 

In the present model, we utilize fixed profiles for 0, O 3 , OH, H0 2 , and 
N0 2 ; these profiles are illustrated in figure 4. The data, for the most part, 
are taken from our large atmospheric simulation above 10 km (recent unpub- 
lished calculations) and are extrapolated to the appropriate ground level 
abundances (e.g., ref. 93). A standard midlatitude ozone profile constructed 
by Krueger and Minzner (ref. 94) has been adopted here. The tropospheric 
water vapor concentrations are derived from a study of midlatitude humidity 
by Sissenwine et al. (ref. 95). Neutral atmospheric densities and tempera- 
tures are taken from the U. S. Standard Atmosphere (ref. 96) for midlatitude 
spring/fall conditions. 

We have performed interactive aerosol-gas calculations with our more com- 
plete aeronomical model and found that, except for the sulfur compounds, none 
of the air constituents is significantly affected by the presence of aerosol 
particles. However, we have not yet treated potentially important radiative 
and chemical feedback effects arising from the optical and catalytic prop- 
erties of the aerosol (e.g., ref. 6), and we have also ignored solution reac- 
tions within the liquid droplets. Some of these problems will be addressed 
in subsequent papers. 
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The lower boundary condition in our photochemical model is specified in 
terms of a fixed boundary concentration, ng^ (molecules cm” 3 ) for each 
species, a, at the ground. A flux necessary to maintain this concentration 
is established at the boundary according to the relation, 

+ = v B (n B - n ) ( 2 ) 

a 

where <j> is the flux (molecules cm ” 2 sec” 1 ); vg is the boundary velocity, 
taken to be 1 cm sec” 1 ; and n^ (molecules cm” 3 ) is the instantaneous species 
abundance. Boundary condition (2) implies that within the model boundary 
layer, counterbalancing surface fluxes seek to maintain n a ~ n^ . The actual 
processes that generate these surface fluxes (emission, absorption, and recom- 
bination) are not specified, but are implied. Except for SO 2 and OCS, whose 
boundary concentrations are taken to be 1.3 ppbv and 0.47 ppbv, respectively, 
each sulfur gas is assumed to have n^ = 0. (For SO 2 , a value of n B of 
1.3 ppbv leads to a surface concentration of 0.6 ppbv in our model.) At the 
upper boundary of the model (which is near the stratopause) , a null flux is 
assumed for each of the sulfur compounds, ignoring a possible small contribu- 
tion of meteoric sulfur to the sulfate layer (ref. 97). Our boundary condi- 
tions are reasonable approximations for a model of the stratospheric aerosol 
layer. 


Predicted Profiles for the Gaseous Sulfur Constituents 

Figure 5 gives our model predictions of typical atmospheric distributions 
of the dominant sulfur gases, including the effects of aerosols. The number 
mixing ratio of OCS is constant in the troposphere — which is indicative of a 
long lifetime in this region — then falls off uniformly above the tropopause 
to maintain an upward flux of molecules and decreases rapidly above 25 km 
where it is photolyzed by sunlight. 

The sulfur dioxide mixing ratio falls off sharply above the ground due to 
washout, which diminishes its abundance to about 0.05 ppbv at the tropopause. 
In the lower stratosphere, SO 2 decreases even more precipitously, being 
trapped and oxidized at these heights. At first, this result might seem to 
contradict the recent observation in reference 55 of a nearly constant SO 2 
mixing fraction for several kilometers above the tropopause level. However, 
these researchers have pointed out that their measurements were made near a 
tropopause break, and they believe that their data were actually taken in 
upper tropospheric air freshly transported into the stratosphere. Since the 
oxidation time of S0 2 in the lower stratosphere is about a week, SO 2 trans- 
ported there would remain intact for some time. Above the aerosol layer, SO 2 
concentrations increase in the region where OCS is photolyzed and again at 
higher altitudes where H 2 SO 4 is photodissociated . 

Sulfuric acid vapor has an interesting distribution in that its lowest 
concentrations occur in the vicinity of the aerosol layer. Obviously the 
aerosol particles absorb H 2 SO 4 , significantly affecting its local abundance. 
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This result suggests that aerosol models that do not include interactive gas- 
particle processes for H2SO4 are probably quite inaccurate in their predic- 
tions for the sulfate layer. 

The H2SO4 mixing fraction increases sharply near the top of the aerosol 
layer (at about 30 km) due to the evaporation of the acid solution droplets 
that rise into this region of warmer, drier air. Notice, however, that the 
total mixing ratio of gaseous sulfur above the aerosol layer (where it is pri- 
marily H2SO4) is smaller by a factor of 2 to 3 than that near the tropopause 
(where it is predominantly OCS) . The gravitational sedimentation of aerosol 
particles effectively acts to reduce the total amount of sulfur in the upper 
stratosphere. 

The sulfur radical HSO3 has a vertical distribution similar to that of 
H2SO4; in fact, it has a larger abundance than H2SO4 over much of the strato- 
sphere. 

The predicted sulfur gas distributions shown in figure 5 are not unique. 
They depend on the assumed model boundary conditions, eddy diffusion profile, 
photochemical rate coefficients, nonsulfurous minor constituent concentra- 
tions, and aerosol abundances; they might also depend to some extent on heter- 
ogeneous physical and chemical processes not yet included in our model. 
Accordingly, the species mixing ratios given in figure 5 should only be inter- 
preted as representative values corresponding to a specific set of aeronomical 
parameters. 


PHYSICAL AEROSOL PROCESSES 


Stratospheric aerosol particles can be formed and existing ones can 
grow in an environment that is saturated with sulfuric acid vapor. Aerosol 
particles can also undergo dynamical motions and interact with each other. In 
this section, we discuss some of the important microphysical processes for 
particles residing in the upper atmosphere, that is, nucleation, condensation, 
evaporation, coagulation, sedimentation, and diffusion, and we describe how 
each process is incorporated into our model. 


Aerosol Droplet Nucleation 

We assume that stratospheric aerosol particles are aqueous sulfuric acid 
solution droplets which are formed in situ by the heterogeneous heteromolecu- 
lar nucleation of H2SO4 and H2O vapors onto condensation nuclei or Aitken par- 
ticles. Reference 98 has shown that the stratospheric nucleation rates for 
homogeneous nucleation or nucleation onto molecular ions are much smaller than 
those for the heterogeneous process. Recent experiments (ref. 99 ) have 
detected undissolved granules in about 30 % of the aerosol droplets collected 
in the stratosphere over a wide geographic range; moreover, most of the other 
droplets, when evaporated, left behind a crystalline matrix suggestive of dis- 
solved core materials. These observations lend credence to the hypothesis 
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that heterogeneous heteromolecular nucleation dominates aerosol formation in 
the ambient stratosphere. Nevertheless, the possibility exists that large 
ions composed of a charged core surrounded by a cluster of polar molecules such 
as H 2 0, HN0 3 , and H 2 S0 4 may be formed in the stratosphere. Some of these ions 
might reach a critical size during their lifetime and grow rapidly into small 
aerosol particles (Castleman, private communication). It is also possible that 
the recombination of large clustered positive and negative ions might lead to 
larger neutral clusters which might act as condensation nuclei. We hope 
to investigate some of these alternative nucleation processes in future work. 

Physically, the heterogeneous heteromolecular nucleation of stratospheric 
Aitken nuclei takes place when sulfuric acid and water vapors can combine to 
"wet" the particle surfaces. Under stratospheric conditions, the newly formed 
particles would consist of the core nucleus surrounded by a (probably impure) 
sulfuric acid-water solution of between 70 and 80% H 2 S 04 by weight, which is 
the range of sulfuric acid concentrations generally found in stratospheric 
aerosol particles (ref. 100 ). 

Because very little is known about heteromolecular heterogeneous nuclea- 
tion, we have parameterized it by assuming that all Aitken particles that are 
mixed into a region where H 2 SO 4 vapor is supersaturated will be nucleated in 
10 6 sec. In NASA TP-1363, we show that this arbitrary choice of nucleation 
time does not significantly affect our model results. 

Stratospheric condensation nuclei are assumed to have a tropospheric 
origin, presumably by atmospheric transport which is simulated in our model 
by eddy diffusion (e.g., see a recent review of stratospheric Aitken nuclei 
in ref. 101). Condensation nuclei may also be deposited above the tropopause 
by volcanic activity, aircraft traffic, and meteoritic ablation, but these 
latter sources have not been incorporated into the ambient model being dis- 
cussed here. 

In determining the concentrations of condensation nuclei in the atmo- 
sphere, we fix their total number density and size distribution at the ground. 
The number density selected is 1200 cm "" 3 and the size distribution varies as 
r _lf for nuclei larger thanO . Olym radius, r being the particle radius 
(0.01 yin is the smallest particle size in the model). The total particle con- 
centration at the ground has been adjusted to give a total of a few hundred 
condensation nuclei and droplets near the tropopause, as has been observed. 

The r -£ * condensation nuclei size distribution is quite different from the 
observed size distribution for stratospheric aerosol droplets (ref. 102). We 
chose the r " 4 distribution because it simultaneously yields the measured 
values of the number of condensation nuclei and their total mass at the 
tropopause level. We have used several other condensation nuclei size spectra 
in our model to determine their influence on the predicted aerosol size dis- 
tribution; our findings are given in NASA TP-1363. 
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Aerosol Growth by Molecular Condensation 


Once sulfuric acid-water solution aerosol droplets are nucleated, they 
can grow by heteromolecular condensation of H2SO4 and H2O vapors. The growth 
of sulfuric acid solution droplets by heteromolecular condensation has been 
discussed in references 34 , 103 , and 104 , to which the reader is referred for 
details. These workers have shown that the growth rate of such a droplet is 
controlled by the rate at which sulfuric acid molecules impinge upon it. The 
fluxes of H2SO4 and H 2 0 vapors to the droplet are balanced so that the droplet 
is always in equilibrium with the water vapor in the environment; that is, the 
H2O vapor pressure of the droplet is always equal to the local partial pres- 
sure of H2O. At a given temperature, the water vapor equilibrium uniquely 
determines the weight percentage of H2SO4 in the solution and, consequently, 
the H2SO4 vapor pressure over the droplet (see ref. 105 ). If the partial 
pressure of H2SO4 in the surrounding air is greater than its vapor pressure, 
some of the sulfuric acid molecules incident on the droplet are absorbed; 
otherwise, sulfuric acid molecules are evaporated from the droplet. The drop- 
let absorbs or evaporates just enough water vapor, essentially instanta- 
neously, to maintain a constant weight percentage of H2SO4 in solution. 


Under these conditions, the growth equation for the droplet is 


dr 

dt 


m D (n - n°) 
a a a a 

rpW(l + X Kn ) 
a a 


X 10 8 


( 3 ) 


oa 


(n - n°) 


where g is the growth rate (yin sec -1 ); m is the mass (g) of an H2SO4 mole- 
cule; r is the droplet radius (pm); p is the droplet solution density 
(g cm -3 ); W is the H2SO4 weight fraction of the solution (we use the terms 
"weight fraction" and "weight percentage" interchangeably; a fraction is 
always implied in our equations); n a is the H2SO4 vapor concentration (mole- 
cules cm" 3 ) near the droplet; and n a is the equilibrium H2SO2 vapor con- 
centration (molecules cm -3 ) over the droplet. (Note that in our calculations, 
vapor pressures have been converted to equivalent molecular concentrations per 
cubic centimeter using the ideal gas law.) The factor of 10 8 in equation ( 3 ) 
is required to convert centimeters to microns. D is a diffusion coefficient 
(cm 2 sec" 1 ) for sulfuric acid molecules in air (ref. 106 ) which depends on the 
mean collision diameter of H2SO4 molecules with air molecules (taken to be 
4 . 3 * 10 _8 cm in this work). The Knudsen number for sulfuric acid vapor diffu- 
sion is given by 

Kn a - t a /r ( 4 ) 

where & a is the H2SO4 molecular free path in air (pm). The parameter X a 
is given by 
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( 5 ) 


1.33 + (0.71/Kn a ) 4(1 - a a > 

X a = 1 + (1/Kn ) + 3ot . 

3. a. 

where a a is the H 2 SO 4 sticking coefficient, which is usually assumed to be 
unity. Note that for large Knudsen numbers, or long molecular free paths, g 
is roughly proportional to a. Hence, the sticking coefficient is an impor- 
tant parameter determining the high-altitude growth rate of solution droplets. 
We discuss the effect of uncertainties in the sticking coefficient more thor- 
oughly in NASA TP-1363. 

Several assumptions are implicit in the derivation of equation (3). The 
effect of dissolved impurities (core material) on the H 2 0 and H 2 SO 4 vapor 
pressures over the droplet is ignored. There are presently very little data 
available to evaluate this effect. In addition, the influence on the tempera- 
ture of the droplet of the latent heat released by vapor condensation is not 
compensated in the growth equation. However, we would expect that strato- 
spheric particles are in close thermal equilibrium with their surroundings. 
Another interaction left essentially untreated is the clustering of water 
vapor to the H 2 SO 4 molecules as they diffuse toward the aerosol droplet. To 
account for this process, an effective cluster diffusion coefficient D a c 
would replace D a in equation (3), and the form of growth equation would be 
modified somewhat (ref. 107). If only one or two water molecules, which are 
small and light, are attached, however, the clustering effect is unimportant. 


According to the preceding analysis of the behavior of H 2 S 04 ~H 2 0 droplets 
in air, we can calculate a fixed set of acid-water solution properties at each 
height in our model in the following manner. At a given altitude, we evaluate 
the equilibrium vapor pressure of water over an H 2 S 04 -H 2 0 solution mixture 
using the empirical relation (ref. 108) 


P°(W) = 


Ci (W) 

(T /T) 1 exp 


(f C e (W)T«-^ 


( 6 ) 


where P£ is the H 2 0 vapor pressure over a flat sulfuric acid solution 
having an H 2 SO 4 weight percentage W; T Q is a reference temperature equal to 
298.15 K; and T is the local air temperature (K) . The constants C(W) are 
tabulated in reference 108 for a discrete set of H 2 SO 4 weight percentages W. 
We selected 30 values of W between 10 and 98% for our model; to obtain vapor 
pressures at intermediate weight percentages, we perform a logarithmic inter- 
polation between adjacent tabulated vapor pressure values. At a given height 
in the model, the equilibrium weight percentage of H 2 SO 4 in solution W is 
determined by matching the calculated H 2 0 vapor pressure P£(W) to the known 
H 2 O partial pressure in the air at the height. Using this value of W, the 
sulfuric acid equilibrium vapor pressure P a (W) is evaluated by an equation 
similar to equation ( 6 ), but with the set of constants C(W) appropriate for 
H 2 S0 4 (also tabulated in ref. 108). For the same value of W, the H^SO^-H^O 
solution density, p (g cm“^) , is obtained from data at 273 K tabulated in ref- 
erence 109. Finally, the vapor pressure of H^SO^ is corrected for the spher- 
icity of the droplets using the Kelvin relation. 
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We have omitted several small effects from the preceding analysis. We 
have not accounted for the Kelvin effect on the H 2 O vapor pressure over small 
droplets. Also, in calculating the Kelvin effect, we have ignored the depen- 
dence of the solution free-surface energy (assumed to be 83 ergs cm -2 ) on the 
solution temperature and composition, and we have not considered the possible 
effects of surface enrichment or volume contraction. Finally, we have not 
extrapolated the acid solution densities below their known 273 K values. 

In the stratosphere, the growth equation gives a negative value for 
dr/dt when the H 2 SO 4 concentration is less than the saturated value n°; 
this corresponds to droplet evaporation. During the evaporation process, 
molecules must escape the free-surface barrier of the droplet, and one might 
logically postulate an "evaporation coefficient" £ which is the counterpart 
of a for condensation. There are no data for a or £ for H 2 SO 4 condensa- 
tion. Nominally, we set a = £ = 1 for aerosol droplets (ref. 110), and 
in NASA TP-1363 we study the effect of these parameters on our model results. 


Aerosol Droplet Cores 

Solution droplets may evaporate to their cores in our model. Accord- 
ingly, we need to distinguish between activated cores — which have a thin 
coating of acid solution and can grow freely by condensation, as can 
droplets — and condensation nuclei, which require a certain nucleation time 
before they can grow (activated cores are equivalent to droplets in our 
model). To do this, we must specify a deactivation time for the activated 
cores which are present in an unsaturated environment, that is, a delay time 
before they become condensation nuclei again. We have arbitrarily assumed a 
zero deactivation time for such cores. As a result, the rate of release of 
condensation nuclei from evaporating particles is only limited by the rate at 
which droplets evaporate to their cores. 

Monitoring the droplet cores in a model is necessary for several reasons. 
First, the size of the cores limits the volume of solution which can evaporate 
from a droplet. Therefore, cores must be considered if one hopes to properly 
conserve the total mass of sulfate in gases and particles. 

Cores also modify the average density of a particle, and the density is a 
factor which influences the rates of particle sedimentation and coagulation. 

If we define f as the fraction of the total droplet volume that is occupied 
by core material, then the average particle density is 

P = fp c + (1 - f)P (7) 

where p is the core density (assumed to be 1.67 gm cm” 3 , as for ammonium 
sulfate). In all our computations of particle settling velocities and coagu- 
lation kernels, we use average droplet densities p corresponding to an 
ambient distribution of aerosol droplet cores which is predicted by our model. 
The average particle densities will change during a perturbation of the 
ambient aerosol layer, but since p Q is close to p anyway, these density 
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changes can usually be ignored (i.e., sedimentation and coagulation rates are 
not significantly altered by small variations in p) . 

To properly evaluate the effect that cores have on aerosol processes, we 
must know the size distribution of the cores within the aerosol droplets. In 
our model we determine the core size distribution using only the first two 
core volume moments. For the cores residing in a collection of droplets of a 
specific size at a given height, we calculate the total core volume and core 
volume squared for those droplets. Our subsequent use of these core volume 
moments to derive the core size distribution is explained in appendix B. 
Although our treatment of droplet cores is an approximate one, we show later 
that it accurately determines core size spectra. 

When a droplet contains several cores — because of coagulation — they 
will probably coalesce to a large degree. Therefore, we would expect to see 
a substantial difference between the size distribution of the cores from which 
a collection of aerosol droplets is formed and the size distribution of the 
cores within the droplets themselves. In NASA TP -1363 we discuss in detail 
the evolution of aerosol core size spectra. 


Water Vapor Growth 

As mentioned earlier, an H2SO4-H2O solution droplet comes into equilib- 
rium with its surroundings by absorbing or evaporating water molecules so 
rapidly that only a negligible amount of H2SO4 is simultaneously transferred 
with the environment. Now imagine a single solution droplet that changes 
altitude. The environmental conditions (temperature, H2O partial pressure, 
etc.) affecting its equilibrium H2SO4 weight percentage and its density also 
change. The droplet absorbs or evaporates water vapor in such a way that 
three quantities are conserved: the number of particles (unity, of course), 

the acid mass of the particle, and the core mass of the particle. Combining 
these three constraints, it can be demonstrated that the particle moves along 
a trajectory in an r-z space, where z is the altitude, such that 

1 _ f o ("r) Jii £n(pw) (8) 

where f Q is the droplet core volume fraction at a reference droplet radius, 
r Q . We refer to the change in the size of an aerosol droplet due to a change 
in its height as water vapor growth. In the stratosphere, the product pW is 
relatively constant over a wide altitude range (e.g., see the data in table 5 
of appendix A). Accordingly, the water vapor growth process is of rather 
limited importance to the stratospheric aerosol. 


dr = _ r_ 
dz 3 


Particle Coagulation 

Coagulation is another important particle growth mechanism which is 
included in our model. Considering a discrete set of particle sizes, the 
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change in the concentration n^ of particles of radius r^ due to collisions 
with particles of radius r^ is (ref. Ill) 


dn. dn- 
x J 

-j— = —r— = -K. .n.n. 
dt dt ij i j 


(9) 


Here, is the coagulation kernel (cm 3 sec" 1 ) for particles of sizes r^ 

and rj. The expression for K-^j which we have adopted is taken from refer- 
ence 111. It might be pointed out that K-. is a function of the diameter 
of the air molecules, which is assumed to be 3.65xl0~ 8 cm, and the kinematic 
viscosity of air, which is evaluated using Sutherland’s equation (ref. 112). 


In our model we consider the coagulation of aerosol droplets with each 
other, and with condensation nuclei, but we ignore self-coagulation of the 
bare condensation cores because of the low probability of their colliding and 
sticking together. A bare nucleus colliding with and adhering to a solution 
droplet is assumed to form a larger droplet. Our aerosol coagulation scheme 
is fully developed in the following two sections and in appendix B. Another 
coagulation process, coalescence, in which larger particles f, run into 11 smaller 
particles due to a difference in their terminal fall velocities, is not impor- 
tant and is neglected in our model. 


Particle Transport by Sedimentation and Diffusion 

The sedimentation of aerosol particles is treated by evaluating the ter- 
minal fall velocity according to the Stokes-Cunningham formula (refs. Ill and 
113). 


The rate of vertical diffusion of particles and gases by atmospheric 
mixing is usually derived from the phenomenological relation 


*D 



( 10 ) 


where cj)£) is the diffusive flux (particles cm“ 2 sec*” 1 ) , D is the eddy diffu- 
sion coefficient (cm 2 sec -1 ) already discussed in the second section of this 
paper. Physical Chemistry of Atmospheric Sulfur, n^ is the air density 
(molecules cm -3 ), and n is the number density of the diffusing species 
(number cm -3 ). Equation (10) states that when particles are uniformly mixed 
in air, their net diffusive flux is zero. We define an "exponentially” 
increasing atmospheric structure parameter 
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°M< ref > 
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where n^(ref) is the air concentration at a reference height. Then equa- 
tion (10) can be written 
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Eddy diffusion is the standard parametric approximation for describing 
the vertical movement of stratospheric gas constituents in one-dimensional 
models. The concept should also apply to particle motions, provided the par- 
ticles are small enough. For large particles, eddy diffusion may be inappro- 
priate; gravitational sedimentation is probably the dominant transport mode in 
this case. Here, all sizes of particles are assumed to diffuse under the 
influence of the same diffusion coefficient, the same one used for gases (see 
fig. 3). We do not expect this to lead to significant errors in the predicted 
transport rates of the aerosol particles. 


Aerosol Removal by Washout and Rainout 

Liquid aerosol droplets can probably be washed from the air very effi- 
ciently by rainfall. We include tropospheric aerosol washout and rainout in 
our model with a removal rate, R^, of about 1 day -1 at the ground, decreasing 
linearly with height (as for the gaseous species) to zero at the tropopause; 
that is, we take 


R w = 


1.2 x 10 5 (1 - z/z T ) sec l 
0 



z > z 


T 


(13) 


where z^ (km) is the tropopause altitude. For a uniformly mixed tropospheric 
aerosol, equation (13) results in an average lifetime against washout and 
rainout of about 2-3 days. We do not include the washout of condensation 
nuclei in our model because we are not attempting to model tropospheric aero- 
sols. In fact, the most important property of condensation nuclei for our 
present purpose is their total concentration at the tropopause. In NASA 
TP-1363 we study the influence of different tropospheric rainout rates on the 
stratospheric aerosol layer. 


Aerosol Parameter Values 

We have tabulated the values of several of the physical aerosol param- 
eters used in our model in appendix A. The vapor pressures, growth rates, 
coagulation kernels, and sedimentation velocities given in appendix A were cal- 
culated using the equations that are referred to in this section and that are 
incorporated into our model. The uncertainties associated with these aerosol 
parameters are briefly reviewed later in table 3. 

In summary, figure 6 depicts, schematically, most of the particle and gas 
mechanisms that are treated in our model. While our simulation is not an 
exhaustive one, we believe that it is sufficiently detailed to give an accu- 
rate representation of the natural aerosol. In the following two sections we 
describe how the aerosol processes delineated above, and shown in figure 6, 
are utilized in our model. 
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THE AEROSOL CONTINUITY EQUATION 


The quantitative effects of chemical and physical processes on aerosol 
particle concentrations can be calculated with a particle continuity equa- 
tion. In this formulation, the rate of change of the particle concentration 
is calculated as the sum of the rates of change due to each individual pro- 
cess. The particle mechanisms which are included in our calculations are 
nucleation, condensation, coagulation, sedimentation, diffusion, and washout. 
Corresponding to these processes, we have the following continuity equation to 
solve: 
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where n is the particle concentration (defined below). 

Various types of aerosol continuity equations have been published before, 
but none of them are directly analogous to equation (14). Accordingly, in 
this section we derive the continuity equations that are required for our 
model. 


Aerosol Variables 

In our model we distinguish between sulfuric acid solution droplets and 
bare condensation nuclei. We define the size distribution of the droplets as 
n(r,z,t), where ndr is the number of droplets per unit volume of air at 
height z and whose radii lie in the range r to r + dr at time t; the 
units of n are particles cm -3 pm -1 . The size distribution of the condensa- 
tion nuclei is defined as n(r,z,t) and has the same interpretation and units 
as n. 


In our model we also need to calculate the first two droplet core volume 
moments; these moments are used to determine the core size spectrum, which in 
turn affects the rates of certain aerosol processes. We denote the total 
volume of the cores in the droplet n(r,z,t) as u(r,z,t) with the units 
ym 3 cm -3 * ym -1 ; u is therefore the summation of n droplet core volumes. 
Because we do not explicitly distinguish between the sizes of the cores in 
individual droplets, we could assume that each droplet has an average size 
core whose fractional volume is 

f (r , z , t) = u(r,z,t)/Vn(r,z,t) (15) 

where V is the volume (ym 3 ) of a particle (either a droplet or a nucleus) of 
radius r. 
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The second core volume moment of the droplets n(r,z,t) is denoted by 
e(r,z,t) and has the units, pm 6 cm“ 3 ym -1 ; e is the summation of the squares 
of n droplet core volumes. In our model e is used to estimate the vari- 
ance, or width, of the core size distribution. 

In the following subsection, we derive the continuity equation for the 
aerosol droplet number density n by developing the corresponding set of rate 
terms for equation (14). Because the continuity equation for n is very much 
like that for n, fewer details of its derivation are given. The continuity 
equations for u and e are obtained by taking the appropriate core volume 
moments of the continuity equation for aerosol droplets having cores of a 
specific size (see appendix B) . Certain terms in the droplet continuity equa- 
tion commute with the core volume moment operators and, as a result, the cor- 
responding rate expressions for n, u, and e all have the same form; these 
"commutative" terms are discussed below. However, the terms for droplet coag- 
ulation, water vapor growth, and evaporation are more involved, and their 
treatment is given in appendix B. For the most part, the continuity equations 
for n, n, u, and e have very similar structures, which simplifies their pre- 
sentation here. 

Many of the quantities that describe physical aerosol processes depend on 
particle size because of weight, drag, and surface tension factors and on 
altitude because of air temperature, density, and trace gas composition 
effects. With this in mind, we do not specifically write out the r and z 
variables in expressions except where clarity demands it. 


The Rate Terms in the Aerosol Continuity Equations 


Heterogeneous heteromolecular nucleation converts bare condensation 
nuclei into solution droplets in situ at a rate which is easily seen to be 


3n\ 

at / 


nuc . 


3n \ 

9t/ 


nuc. 


?), 

*), 


nuc . 


nuc. 


V 

“V 

v v 

V v2 


where is the local nucleation coefficient (sec -1 ). Heteromolecular 

nucleation should not affect the total number of particles (droplets plus 
nuclei), and equations (16) fulfill this requirement. 


(16) 
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The aerosol processes of growth, sedimentation, and diffusion can be 
analyzed by considering the related flux of particles they generate in an 
r-z space; mathematically, the corresponding rate of change of the particle 
concentration due to each process is then the negative divergence of the local 
flux for that process. If we consider aerosol droplet growth by heteromolecu- 
lar condensation, the growth rate, g = dr/dt, from equation (3), can be 
interpreted as a radial "velocity" which produces a radial particle flux, 

4> = gn (17) 

o 

g represents the number of particles per cubic centimeter of air which are 
growing to size r each second. Therefore, the change in <j>g over a small 
radial size interval represents the rate of accumulation or depletion of par- 
ticles in that size interval, or more precisely. 



growth 



(18) 


The growth process also generates radial fluxes of droplet core volume and 
core volume second moment. It is easy to see, however, that cores will be 
affected by growth as follows, for example. 


= - £ (gu) (19) 

growth 

When the sulfuric acid supersaturation ratio n /n ° is less than unity, 
the growth rate g is negative and droplets will evaporate. The extent of 
the evaporation is limited by the cores in the droplets. As we already men- 
tioned, once a droplet evaporates to its core, we assume that it is instan- 
taneously deactivated into a condensation nucleus. Hence, in an evaporative 
environment, there is a net production of condensation nuclei, 

= P (20) 

evap . 

The details of the calculation of P are given in appendix B. 

The evolution of an ensemble of coagulating particles has been studied in 
detail many times (e.g., refs. 26, 111, and 114). The rate of change of the 
concentration of particles of radius r due to coagulation can be expressed 
as the difference between two integrals over the particle size distribution: 
one integral represents the production rate of particles of size r by the 
coagulation of two smaller particles, and the second integral represents the 
loss rate of particles of size r by their coagulation with other particles. 
The coagulation integrals are well known and are not reproduced here. In 
appendix B, however, we do give the coagulation integrals for droplets with 
cores when the core size spectrum is resolved. 
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In our model, aerosol particles are transported by gravitational 
sedimentation and vertical diffusion. The sedimentation flux (particles 
cm” 2 sec -1 pm” 1 ) is 

<P S = -v s n (21) 

where is the settling velocity (cm sec” 1 ); we use the same value of v s 

for all particles (droplets or nuclei) of the same size. For diffusive trans- 
port, the particle flux <f> D is given by equation (12). The particle flux 
equations apply equally to condensation nuclei and acid solution droplets. 

The rate of change of particle concentration due to transport is given 
by the corresponding flux divergence; for example. 
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Often, when the divergence of the sedimentation flux is evaluated in a model 
using standard finite-difference techniques, "numerical" or artificial dif- 
fusion errors can occur. To help eliminate this problem, we have combined 
the diffusion and sedimentation fluxes given by equations (12) and (21) into 
a single composite diffusion flux. For this purpose we define the particle 
scaling factor 
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where D is the diffusion coefficient (cm 2 sec *). Then, it is easy to show 
that 
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and 
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where y is defined by equation (11). 


Since we have assumed that particle transport is independent of particle 
composition, we can readily demonstrate that the droplet core volume flux is 
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(26) 


and that a similar relation holds for s. 
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Droplet water vapor growth occurs with vertical particle displacements. 
Because the vertical diffusion process can, phenomenologically, transport 
cores and droplets at different rates, and because the resulting fractional 
expansion of each droplet is a function of its core size, the water vapor 
growth term is quite complex. Using equation (8), the following expression 
is, however, derived in appendix B: 

®L -*|- * [*—>][•.-■?] I 

growth 


The corresponding water vapor growth equations for u and e are also deduced 
in appendix B. 


Using equations (15) and (24), we can reorganize equation (27) into a 
more useful form; thus, we write 
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with 

g R = f (1 - f >o[^ An(pW)J fcnlyjiCl - f)n] ( 29 ) 

Expressions similar to equations (28) and (29) can be derived for u and e as 
well . 


The rates of tropospheric washout of droplets and their cores are simply 
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Consolidating all of the rate terms we have just discussed, we can write 
out the aerosol droplet and condensation nuclei continuity equations: 
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The continuity equations for u and e have the same general form as equa- 
tion (31). 


The Gas-Particle Interaction 


The growth and evaporation of aerosol droplets can affect the concentra- 
tions of the sulfurous gases in the air surrounding the droplets. In our 
model we specifically account for the interaction of H 2 SO 4 vapor with aerosol 
droplets. The rate at which H 2 SO 4 molecules condense on and evaporate from 
droplet surfaces is determined by the particle growth rate and the total par- 
ticle surface area available. At any altitude, the net rate at which acid 
molecules are condensing onto existing particles is 


dn 
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dt 


lxlO " 12 



47 rr 2 gn dr 


(33) 


where n is the sulfuric acid vapor concentration is Avogadro's number; 
and M is the molecular weight of H 2 SO 4 (the factor 10 ” 12 is needed to con- 
vert units from ym 3 to cm 3 ). Using the droplet growth rate g from equa- 
tion (3), we can rearrange equation (33) into 




(34) 


with 



(35) 


P and L are the sulfuric acid vapor production (molecules cm ” 3 sec” 1 ) and 
loss (sec” 1 ) rates, respectively, due to aerosols. P and L a are included as 
factors in the H 2 SO 4 vapor continuity equation used to calculate its concen- 
tration in air. Our method of solving the continuity equations for the gas- 
eous species has already been outlined in the second section of this paper, 
Physical Chemistry of Atmospheric Sulfur. 


Sulfuric acid vapor condensation on particles, as expressed by equa- 
tions (34) and (35), provides a sensitive feedback mechanism which can mode- 
rate the overall growth rate of stratospheric aerosol droplets. For example, 
if a large number of particles is added to the stratosphere, the H 2 SO 4 vapor 
concentration will be depleted and particle growth will be slowed. On the 
other hand, if the H 2 SO 4 concentration is raised significantly, the existing 
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particles will grow rapidly and decrease the H 2 SO 4 partial pressure. In each 
case, the impact of the perturbation on the aerosol layer is reduced by the 
gas-particle interaction. This interaction is, therefore, an important fea- 
ture of our model. 


A NUMERICAL REPRESENTATION OF THE AEROSOL CONTINUITY EQUATIONS 


Introductory Comments 

We wish to recast the aerosol continuity equations into discrete finite- 
difference relations suitable for computer analysis. The altitude grid 
adopted for the aerosol calculation is the same one used for the photochemical 
analysis; that is, there are 30 equally spaced altitude levels between 0 and 
58 km, with a grid separation h of 2 km. 

Natural aerosol particles usually have diameters spanning several orders 
of magnitude. It is therefore logical, in an approximate model such as ours, 
to utilize a set of geometrically increasing particle size categories such as 
the ones shown in figure 7 (see ref. 26). In this arrangement, particle 
volume doubles from one size category to the next — that is, V 2 = 2Vi , V 3 = 2 V 2 , 
etc. — and the corresponding volume intervals, or bins, associated with each 
size are two-thirds of the category volume. These bin relationships are indi- 
cated in figure 7. For comparison, a "classical" equal-width bin structure is 
also shown in figure 7. Obviously, it takes many more classical bins to cover 
the same size range as a set of geometrical bins. For our model, we have 
selected 25 particle sizes corresponding to particle radii between 0.01 yin 
and 2.56 ym. With the bin boundaries defined in figure 7, the equivalent 
radial interval for each bin is 


r . ( 2 1/3 - l) 2 l/3 

Ar i = ^173 (36) 

where r^ is the radius for the ith size category. Accounting for the bin 
widths, the total aerosol size range that we treat extends from about 0.009 to 
2.9 ym radius. 

In our model, time is divided into discrete steps of t seconds; t is 
automatically increased or decreased during a computation according to the 
maximum rate of change among the gaseous species and aerosol particle concen- 
trations; the same t value is used in both the aerosol and gas phase seg- 
ments of the model. In numerical form, the functional dependence of a varia- 
ble on r, z, or t is replaced by the subscripts i, k, or superscript £ , 
respectively. Our solution scheme only requires that inf ormation be stored for 
two times, t and t + t, or equivalently, £ and £ + 1 ; variables evaluated at 
£ + 1 are, by a previous definition, implicit ones. Our plan here is to 
utilize numerical "implicitness" wherever possible to achieve stability in our 
solutions, while still making the necessary adjustments to satisfy the con- 
servation requirements imposed by the physics. 
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Analogous to the definitions of n, n 
The Aerosol Continuity Equation, we call 


, u, 

n ik’ 


and e 


in the fourth section. 


u £ 

~ ik* ik’ ““ ''ik 

droplet concentration (particles cm 3 ), condensation nuclei concentrations 
(particles cm"^) , droplet core volume (pm 3 cm” 3 ), and core volume second 
moment (pm 6 cm -3 ), respectively, in the radial (or volume) size interval i 
at altitude level k at time £. Hence, for the size-altitude bin (i,k) , 
n ik a s ^ ze “^ nte § rate ^> height-averaged, instantaneous value of n(r,z,t), 

while u^. and are the corresponding (total) core volume moments for the 

droplets in that bin. 


and the solution 


In the following development we limit our discussion, for the most part, 
to the continuity equation for the aerosol droplets, recognizing that there 
is a set of parallel relations for the condensation nuclei and the core vol- 
ume moments. 


The Numerical Treatment of Particle Growth 

To obtain numerical solutions of the aerosol continuity equations, we 
utilize a "discrete” particle formulation for droplet growth and coagulation. 
The discrete approximation assumes that all of the particles in a given size 
bin have the same volume and physical properties; that is, that the particle 
distribution consists of a series of sharp spikes at the bin centers (in vol- 
ume coordinates) . To display the equivalent continuous particle size spec- 
trum, the total number of particles in a bin can be uniformly spread across 
the bin. 

We have also investigated two other numerical approaches. In one, the 
particles ih a bin are assigned a specific size distribution (e.g., uniform 
apportionment across a bin) , and particle growth (and coagulation) rates are 
computed from the integro-dif f erential expressions for these processes by 
insertion of the fixed size distribution. In the second approach, the par- 
ticle size spectrum is assumed to be continuous over the entire particle size 
range of interest, and the radial derivatives and integrals are calculated 
using differencing and quadrature schemes appropriate for a grid of points 
centered in the size bins; in this case, only the size distribution values at 
the grid points are computed. 

For droplet growth, we have modified the discrete formulation using ele- 
ments of the continuous distribution technique — this was done to limit numer- 
ical diffusion in the growth solutions. In the discrete formulation, we 
calculate the rate of volume addition (by condensation) to the particles in a 
certain size bin, and transfer some of these particles to the next largest 
size bin to account for this volume accretion. An accurate estimate of the 
total volume added to particles of all sizes is obtained by evaluating the 
growth integral 
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where 'V is the total particulate volume per unit volume of air (ym 3 cm" 3 ) and 
n-^ jL +1 is a continuous particle size distribution (particles cm"" 3 ym” 1 ). In 
equation (37), we have divided the total particle size range (in our model) 
into radial intervals extending between bin centers; for example, from r^ 
to ^i+i (ignoring the small radial segments at the beginning and the end of 
the size range) . We transfer particles from bin i to bin i-KL according to 
the rate that volume is added in the size interval to r^+j , that is, 

according to the magnitude of dT^/dt. 


We approximate 
law distribution. 


n 


i,i+l 


in each size interval by fitting it with a power 


n. (r) = n. . 

1 , l+l i , l+l 




(38) 


O 

During a time step at a given height, the constants n and a in equation (38) 
are calculated using the initial particle bin concentrations n^ and n^ + ^ 

in the following way; we adjust n and a to simultaneously fit the total par- 
ticle surface areas in bins i and i+1. Thus, we compute two integrals of the 
form, 
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and solve these for n and a. Once n i,i4q( r ) is determined, the rate of 
volume growth between bins i and i+1 is given by the d'V^/dt integral in 
equation (37), and the rate of particle transfer from bin i to bin i+1 by 
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If we adopt a 
size interval 


constant average growth rate, g., for the particles in the 
to r^ » we fi n d from relations (37) to (39) that 
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The use of g introduces a very small numerical dispersion into the particle 
growth solutions. We have imposed a constraint on the value of t. by 
choosing the maximum of the numbers (1) and (n./n-^+i) to evaluate 1 y in 
equation (41). We have done this to enhance tlie computational stability in 
the small particle size range. 
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The net rate of change of the particle concentration in bin i due to 
growth is the difference between the rates at which particles of size i -1 
are entering bin i and particles of size i are leaving. Therefore, to 
complete our development of the discrete aerosol growth equation, we can com- 
bine the (logical) implicit numerical analogues of equation (40) for the size 
bins i -1 and i to obtain 
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There are similar growth equations for the core volume moments u and e. 

As we will show later, our treatment of particle growth effectively 
inhibits the artificial "numerical" diffusion (or rapid growth) of small drop- 
lets into large-size classes in our model solutions. Other workers have also 
addressed this problem (e.g., ref. 115), but their correction schemes are 
often difficult to apply or involve high-order nonimplicit differencing tech- 
niques. In the case of stratospheric aerosols, we always expect to have rela- 
tively smooth and continuous particle size distributions. Hence, the require- 
ments on the growth solutions are less stringent than, for example, in a 
cloud-physics model where very steep size gradients exist, and so our correc- 
tion scheme can be correspondingly simpler. 

According to our previous discussion, the rate of volume addition to all 
of the droplets at height level k at time £ is given by 
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(43) 


In our finite difference solution, equation (43) — which becomes the numerical 
analog of equations (33) and (35) when transformed into H 2 SO 4 molecule units — 
is used to couple the condensation of sulfuric acid vapor on particle sur- 
faces to its gas phase photochemistry. 


Our treatment of droplet evaporation is somewhat different from that of 
droplet growth. In this case, we calculate the time that it would take a 
droplet to evaporate to its average core size, 
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where r c is the average core radius. We use this time constant to estimate 
the rate at which droplets are being converted into activated cores (and 
therefore, condensation nuclei). The conversion rate turns out to be 
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where we have selected an average evaporation rate, | g-^_ 1 | , appropriate for 
the particle size increment to the next smallest size category (one-half of 
the droplet volume would be evaporated in this interval) and where f^. is 
the core volume fraction for the droplets, n^* We ^ ve assumed that, in a 
short time interval, a few M average M droplets in a size bin evaporate all the 
way to their cores — rather than that all of the droplets in the bin evaporate 
to smaller size droplets. Moreover, we have assumed that, on the average, the 
sizes of individual droplet cores do not affect the overall rate of volume 
evaporation. Accordingly, we write the average loss rate of solution volume 
due to evaporation as: 



Equation ( 46 ), which is analogous to equation ( 43 ), is used to calculate the 
net rate of release of H2SO4 molecules during droplet evaporation. It is 
obvious that our evaporation scheme is designed to conserve the number of 
cores and their total volume, as well as the mass of sulfuric acid, for a col- 
lection of evaporating droplets. 


A more detailed evaporation scheme would allow droplets to evaporate into 
smaller size bins, carrying along the smaller cores and releasing the largest 
ones as condensation nuclei. Under stratospheric conditions, however, any 
acid solution droplets that are transported into a region where evaporation 
occurs will usually disappear very rapidly. In addition, since we do not know 
the exact size distribution of the cores within the droplets, it is difficult 
to determine accurately the fraction of the droplets that would release cores 
during an incremental decrease in the droplet size. Because a detailed^evapo- 
ration treatment is too involved, we utilize the core volume moments, u^j 
and e^, to estimate the core size distribution and use this approximate 
distribution as a guide in placing evaporated cores into appropriate condensa- 
tion nuclei size bins (see appendix B) . In the numerical formulation of 
evaporation, as outlined in appendices B and C, we calculate the rates of 
production of various sizes of condensation nuclei from evaporating drop- 
lets — the p|k for equation ( 32 ) — using the droplet evaporation rate from 
equation ( 45 ) and the core size distribution described in appendix B. 


Since the vapor pressure of a droplet increases with decreasing radius 
(Kelvin effect), an environment which is only marginally supersaturated with 
H2SO4 may have the large particles growing while the small ones are evaporat- 
ing. Under these conditions, special care must be taken in handling the par- 
ticle growth rates for the size bins where the transition from growth to 
evaporation occurs. 
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The Numerical Solution for Particle Coagulation 


Coagulation must be treated approximately in a numerical model because of 
the complexity of the integrals describing this process (see the discussion in 
appendix B) . At the beginning of the subsection entitled The Numerical Treat- 
ment of Particle Growth, we mentioned three possible techniques which might be 
used to calculate particle coagulation rates: one technique is based on a 
discrete particle size spectrum and the other two employ continuous particle 
size distributions. For coagulation, these three numerical schemes lead to 
finite-difference equations having the same general form. Even so, the two 
techniques which assume continuous size distributions can have problems 
related to particle number and volume conservation. Actually, all three solu- 
tions give about the same results for the average rates of particle coagula- 
tion. Since the discrete formulation is much simpler to use, however, we have 
adopted it here; later we demonstrate that the discrete technique is also 
quite precise. 


For each pair (i,j) of particle size bins in our model, a coagulation 
kernel is calculated at each height using the equations referred to in the 
subsection entitled Water Vapor Growth: 


K. M = K(r . ,r . ,z ) = K 
ljk i J k jik 


(47) 


The coagulation kernel is symmetric in the size indices i and j . The data 
tabulated in appendix A indicate that the coagulation kernels can vary up to 
40% between adjacent size bins. Therefore, our numerical solutions of coagu- 
lation effects may contain some additional numerical dispersion due to our 
use of these average kernels. 


In our discrete numerical formulation, we consider coagulation to be an 
interaction process between pairs of particles whose sizes are limited to a 
finite set of discrete values. Thus, instead of analyzing the continuous 
coagulation integrals, we proceed directly to a discrete analog of those inte- 
grals. Kritz (ref. 26) was one of the first workers to adopt geometrically 
increasing size bins to study discrete particle coagulation; this approach 
leads to some interesting numerical interpretations of the coagulation pro- 
cess. First consider the loss rate of particles of size i by coagulation 
with all other size droplets. Clearly, whenever a particle i collides with 
and adheres to a larger particle, particle i disappears from its own bin; 
the rate of disappearance is (suppressing the time and altitude variables) : 


^coag. 
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For encounters with smaller particles, however, the rate of removal of one i 
droplet from a size bin would require several sticking collisions with the 
smaller particles to add enough volume to make a particle of the next larger 
size; in fact, it would require Vj^/Vj collisions. This is an interpretation 
of coagulation as a discrete growth mechanism. On the other hand, if we 
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consider the intermediate size coagulated particle composed of droplets 
i and j, we can subdivide it into a fraction Vj/V^ of an i+1 droplet and a 
fraction (1 - Vj/V^) of an i droplet. That is, a coagulation event which 
increases the volume of an i droplet by less than a discrete size increment 
can be interpreted as transferring only a fraction of the i droplet to the 
next size bin, while preserving particle number and volume in the process. In 
either case, the net rate of loss of droplets i by coagulation with smaller 
droplets is 
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For self-encounters, the kernel K. . is not multiplied by a "symmetry" factor 
since the effect of counting self-collisions twice is already compensated in 
equation (49). Combining equations (48) and (49), the total loss rate of 
droplets i by coagulation is 

= -n . L . (coag . ) 

coag . 
loss 

+ E 

Lj<i j>i 


K. .n. 
iJ J 


(50) 



Before proceeding, we point out that the coagulation terms all contain 
products of aerosol concentrations; we cannot easily solve the nonlinear 
finite-difference equations that result when each of the concentrations is 
evaluated implicitly at the time, £+1. Accordingly, we make the concentration 
products "implicitly" linear by evaluating n^ at £+1 and nj at £ . 

(This mixing of time dependencies actually results in a more precise numerical 
solution than one which evaluates the coagulation rate at the beginning or the 
end of a computational step.) In finite difference form, the coagulation loss 
rate in equation (50) becomes, 
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The rate of production of particles i+1 by coagulation must equal the 
rate of loss of particles i in collisions with particles j<i; conversely, 
it must also equal the sum of the rates of loss of particles j<i in colli- 
sions with particles i weighted by the volume factors Vj/V^ to account for 
the actual number of collisions which would be required to form a particle 
i+1. As before, an equivalent argument is based on the formation of an inter- 
mediate size coagulated droplet i+j of which a fraction V j/ V i is assigned 
to bin i+1. Following our convention for implicitness, ,! we can numerically 
calculate the coagulation production rate of particles i+1 as an average of 
the coagulation loss rates for particles i and for particles j<i, namely. 
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The 0j_j factor is necessary to eliminate double summing when self- 
coagulation terms are evaluated. Equation (52) is balanced "implicitly" and, 
in conjunction with equation (51), conserves discrete droplet volume exactly. 


We can include the effects of coagulation between condensation nuclei and 
aerosol droplets in our formulation once we make the following assumptions: 

(1) that nuclei are lost and droplets are formed during the coagulation pro- 
cess; (2) that the nucleus-droplet coagulation kernels are identical to those 
for droplet-droplet collisions (ignoring differences in the sticking effi- 
ciencies and particle densities); and (3) that each nucleus contributes its 
entire volume to the resulting droplet core. With these guidelines, we can 
derive the relations: 




3t ) 


coag . 
loss 


8n £ 

ik 


3t /. 


) 


coag . 
loss 


* £+1 
*ik 


^ , K ijk n jk 

{all J} 


(54) 



K ijk vT (“5k + “ 


jk 


) + X 


l 

l jk 


j>i 



(55) 


34 



( 56 ) 


9n 


ik 


9 1 


a * 

9n.. i 
l+l ,k 

9 1 


coag. 
prod . 


coag . 
prod . 


,a+i __ (, v j\ i 

n ik ijk \ v ± / n J k 




1 V' „ / V j \ / l+l * . t 

2 / AM/rj v’ik n jk + "ik n 


£+1 £+1 ^ & 

jk + n ik n jk 




, & ~ £+1 ~ £+1 £ 

+ n n .. + n .. n + 

ik jk ik jk 


£+l\ a 

n ... n ... 10.. 

ik j k / ij 


(57) 


Particle coagulation at any altitude must conserve the total particle 
volume. To check this in our equations, we compute, at altitude level k. 
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Accounting for the possible coagulation of particles out of the largest size 
bin, we can show, after some algebra, that the sum on the right is, as 
required, identically zero. 

The effect of coagulation on the evolution of droplet core volume moments 
is discussed in appendix B. The resulting coagulation equations have the same 
general structure as those for the droplets. The radial boundary conditions 
for droplet growth, evaporation, and coagulation are described in appendix C. 


The Numerical Treatment of Particle Transport 


The sedimentation and diffusion terms in the aerosol continuity equation 
are handled as follows. The composite particle flux from equation (24) is 
numerically evaluated midway between altitude levels, for example, as 
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Then we calculate the flux divergence. 



ocu. . 

diff . 


+ 


1 r SD _ SD 

h [_^±,k+(l/2) <P i,k-(l/2) 


(60) 


35 



In equation (59), the expression for D/yy is a logarithmic interpolation to 
the midaltitude point. One can see by inspection that our transport formula- 
tion for sedimentation and diffusion conserves particle mass exactly in a 
vertical column (accounting for the boundary fluxes) . 


Some Additional Numerics 

Water vapor growth, expressed by equations (28) and (29), can be analyzed 
in the same way as growth due to acid vapor condensation. If we calculate the 
rate of volume addition to the droplets in size bin i due to water vapor 
growth, using the assumption of a discrete particle size spectrum, and if we 
adopt a relation similar to equation (40) to describe the growth process, we 
can estimate the rate of particle transfer from bin i to bin i +1 using 
the time constant 


"5F ' [ B 37 

T ik 

where y, defined in equation (41), limits numerical diffusion effects. The 
finite difference equivalent of equation (61), which we use in our model, 
except at the boundaries, is 

(l - f 5 ' 

i , k+1 \ i,k+l 

(l - f 5, 

i,k-l\ i,k-l 

The corresponding water vapor growth rate terms for u and e are very similar 
to equation (62). 

In connection with aerosol growth processes and our numerical analogs of 
them, we note that we have not written a separate continuity equation to keep 
track of the total mass of sulfuric acid within the aerosol solution droplets. 
Rather, we have carefully conserved the total droplet volume and core volume 
(and, hence, the acid solution volume) during each aerosol interaction. For 
proper H 2 SO 4 conservation, this requires the inclusion of water vapor growth 
effects in our solutions. Because the water vapor growth effect is always 
small, however, we only need an approximate numerical representation for it — 
one such as equation (62) affords. In our model, we have also ignored water 
vapor growth in the upper and lower boundary layers (k = 1 and 30) and in 
regions where droplets are evaporating. The latter omission is justifiable 
because according to our model, aerosol droplets do not undergo large vertical 
displacements as they are evaporating (the evaporation rates above the aerosol 
layer are quite large) . 

The heterogeneous nucleation and washout rates in equations (16) and 
(30), respectively, can, quite simply, be replaced by their appropriate aver- 
age values for each size-altitude bin; for example, R^n ^ 1 . The particle 
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concentration in each case is evaluated implicitly at time £+ 1 , which happens 
to couple the condensation nuclei and aerosol droplet continuity equations; a 
technique to solve the coupled equations simultaneously is given in appen- 
dix C. The heterogeneous nucleation rate R^. is lxl0 ~ 6 sec -1 whenever the 
sulfuric acid supersaturation ratio n a /n a ° is greater than 1 in the model; 
otherwise it is zero. The droplet washout rate is specified in equation (13). 


Recapitulating our previous development, we can write the continuity 
equation for an evolving aerosol system in finite-difference form as, 
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where we have replaced 3n/3t on the left-hand side of equation (14) by its 
numerical equivalent. There are similar equations for the condensation nuclei 
concentration and the core volume moments. Some details of a computer algo- 
rithm to solve these equations are given in appendix C. 


MODEL VALIDATIONS AND SIMULATIONS 


In this section we discuss, first, the numerical validation tests used to 
evaluate the accuracy and behavior of our model finite difference algorithms 
and, second, the ambient stratospheric aerosol simulations made for compari- 
son to observational data. We have also summarized in tables 2 through 4 all 
of the major assumptions involved in the construction of our model. 


Numerical Validation Tests 

To validate the physical aerosol mechanisms as they are formulated in our 
model, several simple numerical calculations have been made. In each case, 
all aerosol processes except one are shut off and, under these controlled con- 
ditions, the resulting model prediction is compared to an "exact” analytical 
solution. 
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To obtain useful analytical solutions of the aerosol continuity equa- 
tions (31) and (32), certain simplifying assumptions can be made. For 
example, at a fixed altitude in the sulfate layer we can consider the steady- 
state continuity equation (31) in the form 

^t[g(r)n(r)] + = P(r) (64) 

R 


where x^(sec) is the local aerosol residence time — including the effects of 
sedimentation, diffusion, and coagulation — and P (droplets cm -3 sec” 1 ) is 
the local aerosol production rate, due to nucleation and coagulation, with 
possible contributions from other sources such as aircraft engines. The solu- 
tion of equation (64) is 
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An especially simple solution of equation (65) results when we set P = 0 
and ignore the dependence of g and t r on particle radius; that is. 


n (r ) = n(r Q )exp[-(r - r^/r^] (66) 

where 

r m " gT R <67 > 

is the characteristic radius for this exponential size distribution. The pre- 
dicted average particle size, in this case, is determined by the average 
radial growth that occurs during the droplet f s residence time in the upper 
atmosphere . 

A numerical aerosol growth experiment is illustrated in figure 8. In 
this calculation the number of particles in the first size bin is fixed 
(1 particle), a constant growth rate g of 0.5 pm yr” 1 is specified, and all 
other aerosol interactions are deactivated except for an assumed removal pro- 
cess with an exponential lifetime t r of 0.2 yr . The corresponding numerical 
and theoretical steady-state particle distributions for these conditions are 
presented as curves at time t = 0 in figure 8; the analytical solution is 
represented by equations (66) and (67). The difference between the two solu- 
tions can be roughly interpreted as a 5% underestimate of the particle growth 
rate in the numerical simulation. This result should be compared to another 
model prediction obtained without a correction for numerical diffusion in the 
particle growth rates. In this instance, the large (~~2.5 pm) particle con- 
centrations were overestimated by several orders of magnitude. Numerical dif- 
fusion normally occurs because at each time step some of the particles in each 
size bin are allowed to grow across the entire bin. We have suppressed 
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numerical diffusion by limiting the number of these "fast-growing" 
particles — by carefully calculating the rate at which volume is being added 
to each bin (see A Numerical Representation of the Aerosol Continuity Equa- 
tions) . 

When the droplet growth rate is doubled in our growth experiment, the 
ensuing development of the particle size distribution is that shown in 
figure 8. Our model solution for the time-dependent problem requires an aver- 
age of only about 30 time steps per year. Note that the agreement between the 
analytical and numerical solutions is quite good, except for the very small 
systematic underestimate of the growth rate. 

The results of the growth experiment clearly demonstrate that "numerical" 
diffusion is not polluting the large particle size range with artificially 
grown droplets in our model. We have also checked to see that the number of 
particles is conserved exactly during the growth process. 

Figures 9 and 10 show the results of a particle coagulation experiment 
in which all of the other aerosol processes have been deactivated. Starting 
with one particle in the first size bin, we follow its coagulation history 
with our numerical coagulation scheme. Time is expressed in units of the 
characteristic coagulation time, t q (e.g., see ref. 116) and the experiment 
is continued for 10 time units. The coagulation time scales of importance to 
stratospheric aerosol problems are usually much less than 10 t q , since t q 
can be as long as several months in the lower stratosphere. For our calcula- 
tion, only 20 equi-spaced time steps are taken for each unit of t q . By com- 
parison to the exact solution of Smoluchowski (e.g., ref. 116), the model 
predictions in figure 9 are excellent. The data in figure 10 illustrate the 
corresponding development of the coagulation size spectrum and show that it 
is also in very good agreement with the exact Smoluchowski spectrum even after 
a long relaxation period. (The particle distributions in figure 10 are not 
size distributions, but are bin-population distributions for our geometric 
size bins. Accordingly, one of our size bins may contain a number of Smolu- 
chowski bins — see fig. 7 — a fact we have taken into account in computing 
the equivalent Smoluchowski populations.) 

To test our numerical formulation of sedimentation and diffusion, we have 
made the calculation shown in figure 11. In this experiment the atmosphere is 
assumed to have a constant scale height of 10 km, and both the settling veloc- 
ity and diffusion coefficients are made inversely proportional to the air con- 
centration. A point source of particles is placed at 20 km, and the aerosol 
concentration at 12 km is forced to be a very small value. All other particle 
interactions are turned off. 

Analytical solutions for this problem can be obtained using a simplified 
version of the particle continuity equation (31) with only sedimentation and 
diffusion of a monodispersed aerosol, but with an added source of particles, 

Q (particles cm ~ 3 sec" 1 ). The solution of the resulting steady-state equa- 
tion gives the particle concentration profile. 
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( 68 ) 
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<(>p is the net (downward) particle flux (particles cm 2 sec 1 ) which is 
related to the particle source Q by 
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For a point source of particles at a height z p the particle mixing ratio 
ny decreases exponentially above z p with a scale height D/v s where v g 
is the settling velocity. The analytical solutions for our settling experi- 
ment are plotted in figure 11 along with the model predictions for two dif- 
ferent diffusion rates and two assumed particle sizes. Obviously, the 
numerical model gives a very accurate representation of one-dimensional aero- 
sol transport processes. 

Note in figure 11 that the vertical profile of the large (1.0 pm) parti- 
cles is very sensitive to the vertical diffusion rate. Large particles might 
therefore be a useful tracer of vertical atmospheric motions. On the other 
hand, the steep vertical gradients in the large particle concentration pro- 
files can lead to numerical diffusion errors in their solutions. In another 
experiment, where sedimentation and diffusion fluxes were calculated sepa- 
rately (i.e., without utilizing the combined flux of equation (24)), the pre- 
dicted 1.0-ym particle concentrations at 30 km were overestimated by more 
than an order of magnitude. Obviously, our use of a combined particle flux 
(defined by eq . (24)) results in an overall diffusion-like transport process 
in which numerical diffusion effects are minimal. 


The vertical dispersion of a narrow layer of aerosol particles due to 
the eddy diffusion coefficient in figure 3 is illustrated in figure 12. 
Starting with a concentration of the 100 particles cm“ 3 in the 2-km-thick 
altitude interval centered at 20 km, the particles are transported to other 
heights only by diffusion. We have verified that the diffusion process 
exactly conserves the total number of particles. Notice that the strato- 
spheric residence time of the aerosol due to diffusion is about 2 to 3 years 
in this model, although it is obviously a function of the assumed eddy diffu- 
sion profile. 

We have also checked the accuracy of the numerical algorithms in our 
model in other ways: by confirming that the coagulation process exactly con- 

serves core mass and particle volume, that the overall atmospheric budget of 
sulfur atoms in gases and particles and of core material in nuclei and drop- 
lets correctly balances, and that the total number of nuclei and droplets is 
properly accounted for. 
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Numerical Model Simulations 


We have made steady-state model predictions of the properties of the 
ambient stratospheric aerosol layer; these predictions of its characteristics 
are compared to experimental data in a series of figures (figs, 13 to 17). To 
make our calculations, we have used the aerosol model exactly as it is 
described in the preceding text. The values of the physical parameters and 
the techniques of numerical solution are those given in earlier sections of 
this paper and in the appendices. 

Interestingly, even though we have not attempted to "tune* 1 our model by 
varying any of the "adjustable 11 physical parameters — and there are many of 
them — the predictions for the five aerosol characteristics shown in fig- 
ures 13 to 17 are reasonably good. We should also point out that there are 
many uncertainties in our calculations, as the comments in tables 2 through 
4 indicate. Moreover, natural atmospheric variability and inhomogeneity make 
it difficult to compare steady-state, one-dimensional model predictions with 
a limited number of observational points. Nevertheless, the overall structure 
of the stratospheric aerosol layer is well established, and a one-dimensional 
model can be a useful diagnostic tool for studying the general properties of 
the layer. In NASA TP-1363, we present a detailed analysis of our model pre- 
dictions in relation to the observed characteristics of the aerosol layer and 
the known physical processes of aerosol particles. In NASA TP-1363, we also 
give the results of a comprehensive study of the sensitivity of our model pre- 
dictions to the uncertainties in the model parameters. Therefore, we will 
conclude this paper with only a brief description of the model calculations 
shown in figures 13 to 17. 

The sulfate mass mixing ratios plotted in figure 13, for example, indi- 
cate that the model gives roughly the same amount of sulfate mass in the aero- 
sol layer as has been observed. This result suggests that the stratospheric 
sulfur balance in our aerosol model is probably close to reality. 

Predicted and measured total particle mixing ratios in the stratosphere 
are shown in figure 14. The scatter in the experimental data could be indica- 
tive of atmospheric variability, which might be partly related to the rate of 
vertical particle transport in this case; for Junge’s data, the scatter is 
probably experimental. The steady-state model predictions presented here are 
believed to represent an average state for the ambient aerosol layer. How- 
ever, we plan to study the effects of variability on the observable charac- 
teristics of the layer with our time-dependent model at a later date. 

Figure 15 illustrates that our calculated concentrations of large parti- 
cles (radii > 0.15 pm) are in accord with a number of stratospheric in situ 
measurements. Figure 16 shows that the concentration ratios of the large par- 
ticles (those with radii > 0.15 pm to those with radii > 0.25 pm) are also in 
agreement with observations. In both cases, there is a substantial range of 
measured values, with the model calculation being generally representative of 
the data. 
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The computed stratospheric particle size distributions at 16 and 20 km 
are compared to several empirical size distributions in figure 17. The model 
reproduces the observed stratospheric aerosol zold size spectrum quite well. 

To validate an aerosol model by comparing its predictions with atmo- 
spheric observations requires a diversified and high-quality data base. 
Unfortunately, existing aerosol data are not completely satisfactory for this 
purpose. In NASA TP-1363, therefore, we utilize our model to explore the pos- 
sible experimental studies which might be most helpful to atmospheric scien- 
tists attempting to understand the origin and composition of the aerosol layer. 


CONCLUSIONS 


We have developed a new time-dependent, one-dimensional model of the ter- 
restrial sulfate aerosol layer. Our model treats atmospheric photochemistry 
and aerosol physics in detail and includes the interaction between gases and 
particles explicitly. We have demonstrated that the numerical algorithms in 
our model are quite precise and have shown that the simulated aerosol physics 
generates a particle layer with most of the observed properties. In NASA 
TP-1363, we discuss at length the sensitivity of the calculated sulfate layer 
properties to changes in a large number of aeronomic and aerosol parameters. 

Our model can predict almost all of the observable aerosol characteris- 
tics, including the height and shape of the layer, the total mass and number 
mixing ratios, the large particle mixing ratios, the small-to-large particle 
ratios, the acid solution droplet composition and size distributions, and the 
droplet core volumes. Moreover, it can simultaneously predict several mea- 
surable aeronomical quantities, such as the abundances of gaseous sulfur com- 
pounds in air, the vertical profiles of neutral trace constituents like 
methane and nitrous oxide, and the stratospheric residence times of gases and 
particles. As a result, the model calculations are highly constrained by 
data and are therefore useful for experimentation in parameter sensitivity. 

Our aerosol model may be used in any one of several different computa- 
tional modes: as a steady-state model to study average long-term aerosol 

characteristics or to perform aerosol sensitivity tests; as a time-dependent 
model to investigate perturbations of the sulfate layer (man-made or natural), 
short-time scale variability of the layer, or aerosol aging; or as a limited 
physical model, after some simple engineering alterations, to look at indi- 
vidual aerosol processes such as coagulation or growth, or to estimate the 
removal rates of particulate matter from the stratosphere. Our aerosol simu- 
lation has one severe limitation — it does not account for the horizontal 
spreading of gases and particles. We plan to correct this deficiency by 
developing a two-dimensional version of the present model based upon an exist- 
ing two-dimensional model of stratospheric trace gases (e.g., see refs. 120 
and 121). In the meantime we expect to apply our one-dimensional model to 
some unresolved aerosol-related problems including: the identification of the 

mechanisms that form stratospheric aerosol particles from gases and nuclei and 
which determine their eventual size spectrum; a clarification of the role of 
volcanoes and other sulfur sources in maintaining the sulfate layer; and an 


42 



estimation of the effect of man-made emissions of sulfur gases and solid 
particulates on the ambient aerosol layer. 


Ames Research Center 

National Aeronautics and Space Administration 

Moffett Field, California 94035, November 29, 1978 
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APPENDIX A 


PHYSICAL AEROSOL PARAMETERS 


In this appendix we have tabulated values for several of the aerosol 
parameters discussed in the third section of this report. Physical Aerosol 
Processes, For sulfuric acid-water solution droplets of about 0.3 pm radius, 
table 5 lists the solution density, the equilibrium H 2 S0 4 vapor concentration 
over the solution, and the weight percentage of H^SO^ in the solution at the 
altitude levels used in our model, up to a height of 40 km. The atmospheric 
temperatures in table 5 are taken from reference 96 for midlatitude spring/ 
fall conditions and the water vapor concentrations from figure 4. Table 6 
gives the values of the aerosol growth "kernel" g oa (defined by eq. (3)) for 
the entire range of particle sizes used in our model and for altitudes up to 
40 km. For brevity, only every third particle size has been selected for 
table 6 (and tables 7 and 8 as well). Droplet coagulation kernels at a height 
of 16 km are listed in table 7, and aerosol settling velocities are presented 
in table 8. 
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APPENDIX B 


AEROSOL DROPLET CORES: THEIR SIZE DISTRIBUTION AND EFFECT ON AEROSOL BEHAVIOR 


If we consider a collection of aerosol droplets that have identifiable 
cores, we can categorize the droplets by their total volume v(ym 3 ) and 
their core volume w(ym 3 ), where w < v, Likewise, we can define the droplet 
size distribution q(v,w,t), where n dv dw is the concentration of droplets 
per cubic centimeter of air whose volumes v lie in the range, v to v + dv, 
and whose core sizes w lie in the interval, w to w + dw, at time t; the 
units of n are particles cm™ 3 ym" 3 ym -3 . (Note that the spatial dependence 
of n is suppressed and that particle volume rather than particle radius is 
used as a size coordinate. Because the transformation between volume and 
radius is straightforward, v and r are easily interchangeable size parameters 
in our equations.) 


The aerosol variables employed in the text can be obtained from q by 
integration over core sizes: 
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n(v,t) = I dw n(v,w,t) 
'o 
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u(v,t) =1 dw w q(v,w,t) 

*o 


(Bl) 


(B2) 


e(v,t) = I dw w 2 q (v,w, t) (B3) 

"o 

These equations define the total core volume "moments" u and e. 


We can readily formulate a continuity equation for the droplet distribu- 
tion n; it will contain a rate term(s) for each physical aerosol process in 
our model. By taking the appropriate core volume moments of the q conti- 
nuity equation, we can derive the continuity equations for n, u, and e. The 
solutions of these equations are discussed at length in the text. Calculated 
values of u and e are used to roughly characterize the droplet core size 
distribution; the manner in which this is done will be described shortly. 


In the q continuity equation, we can show that the mathematical opera- 
tors associated with particle growth, sedimentation, diffusion and washout all 
commute with the core volume moment integrals. As a result, the equivalent 
continuity equation terms for q, n, u, and e are all similar, anc the dis- 
cussion of the n terms in the text can readily be generalized. However, we 
do have to consider more carefully the aerosol processes which are core depen- 
dent: that is, coagulation, which redistributes cores; water vapor g' owth, 

which depends on core fractions; and droplet evaporation, which is limited by 
core size. 
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DROPLET COAGULATION 


As a starting point for our analysis of aerosol coagulation, we write the 
coagulation equation for droplets with cores. 




^3 1 * ^ = Y f dv'K(v' ,v - v') / dw'n(v' ,w' ,t)n(v - v' ,w - w',t) 

* r\ 


- r)(v,w,t) f dv'K(v,v') f dw' n (v T ,w' , t) 

n " n 


(B4) 


where K is the coagulation kernel. The limits on the dw 1 integration in 
the first term of equation (B4) must, in general, be modified to satisfy the 
following constraints on the core volumes: 


w T < w \ 

w - w T < v - v T J 


(B5) 


Equation (B4) is a generalization of the well-known integral equation describ- 
ing particle coagulation. The 1/2 factor in the first term of equation (B4) 
arises because the integration over the entire range of dv' is found to 
count each pair of particle collisions exactly twice when one carries out the 
dw f integral indicated. 

Taking the first three core volume moments (i.e., 1, w, w 2 ) of equa- 
tion (B4), the time derivative and coagulation loss integral reduce to: 
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The coagulation production term, however, contains the core integrals, 
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with m = 0,1,2. Transforming the iterated integral (B7) to the variables 
jw r ;w" = w - w T | and using the fact that the physical limits on w* and w" 
are 0 to v 1 and v - v T , respectively, expression (B7) becomes 
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In this form, the core integrations can be factored, and the resulting 
coagulation production integrals for the core volume moments of interest are: 

v f ,t) 

v r ,t) + u (v - v f , t)n(v 1 , t) 
v T , t ) + 2u(v T , t)u(v - v T , t ) 

+ e(v - v , ,t)n(v , ,t) (B9) 

The form of expressions (B9) indicates that we have to treat the coagulation 
terms for u and e with special care in our numerical formulation (the coagu- 
lation terms for n are already developed in the text). Interestingly, the 
core volume moment equations are not coupled to higher order moments, and we 
can solve them without recourse to a closure approximation. 

In a finite difference scheme which is based on a discrete set of droplet 
sizes, there is some ambiguity about the fate of the cores in coagulating 
droplets. Since our particle volumes double between size bins, several small 
droplets must, in general, coalesce to produce a droplet of the next largest 
discrete size; the resulting core is then composed of all the individual coag- 
ulated cores. However, a simpler, and more logical interpretation of core 
coagulation — which is also consistent with our discrete numerical treatment 
of droplet coagulation — takes into account only the rates at which pairs of 
particles are coagulating. In most cases, droplets of sizes i and j < i 
coalesce into a particle whose volume + Vj is intermediate between the 

discrete sizes and V-^ +1 . The average core volume (denoted here by W) in 

this intermediate droplet is given by 
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dv' K(v ' , v - v 


{ n(v' , t)n 
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In our droplet coagulation scheme, we divide such an intermediate size 
coagulated droplet into two fractional droplets of sizes i and i+1 (while 
conserving the total droplet volume and number — see the subsection. The 
Numerical Solution for Particle Coagulation, in the main body of the text). 
However, we now want to also specify that the average core volume fraction f 
for each of the droplet fragments should be the same as for the intermediate 
droplet. That is, we require 


V. + V. 
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W. 
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V. 
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w i + l 

v ifl 


(Bll) 


where W. and W i+1 are the volumes of the cores assigned to the fractional 
droplets. The requirement of a constant core volume fraction during coagula- 
tion is a constraint that can be justified by the following thought experi- 
ment. Assume that we have an initial ensemble of droplets, all with the same 
core fraction. Then, as the ensemble evolves by coagulation, every droplet 
must always have that same core fraction. Equation (Bll) automatically satis- 
fies this condition. 


il 
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The rate of coagulation of the droplets i and j is Ky n^iij , but only a 
fraction, Yij - Vj/V^, of a large droplet (i+1) will be formed by each coagu- 
lation. Accordingly, the rate of core volume transfer to the droplets of size 
i+1 is the product of W i+1 , K^nyuj and Yij » or 


p i+i < u > 


coag. 




+ u .n . )0 . . 
3 i iJ 


(B12) 


where 0=1 except when j = i, in which case 0=1/2 (again, see The 
Numerical Solution for Particle Coagulation). Equation (B12) is the discrete 
analog of the core volume integral (B9) . During coagulation, a fraction, 

1 - » of an i droplet is also formed. This partly offsets the loss of 

the original i droplet which coagulated, but its core has been modified 
somewhat. The net core volume change in bin i is 



(B13) 


Equation (B13) shows that, because of core rearrangement, the coagulation of 
droplets i and j leads to the transfer of core volume from droplets j to 
droplets i at the rate 


AP. (u) 
i 


coag. 




(B14) 


Moreover, core rearrangement also alters the net rate at which droplets i 
lose core volume by direct coagulation; without rearrangement, the direct loss 
rate (the u^ term in eq. (B13)) would be Kj_j Yj_j u i n j • For the smaller drop- 
lets j, however, the direct core loss rate (eq. J (B6) ) is unchanged and there 
is no core rearrangement. By summing the core volume production and loss 
terms for the droplets j, i, and i+1, it can be shown that the total core 
volume will be conserved during coagulation. 

Our coagulation scheme attempts to follow the evolution of the cores 
within a collection of droplets of discrete size. Because of the finite width 
of the droplet size bins, it is reasonable to assume that collisions between 
the particles in a small bin and those in a large bin will affect the distri- 
bution of the cores in the larger bin; that is, that some of the collisions 
will produce droplets with a modified core structure within the boundaries of 
the larger bin. Our numerical approach to droplet coagulation, as outlined 
above, is actually a discrete representation of the exact coagulation integral 
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(B9) — one that is limited, however, by the lack of a clear definition of the 
droplet size distribution in each finite size bin. 


The effect of coagulation on the second core volume moment e is calcu- 
lated in a manner similar to that used to calculate the first core volume 
moment u; that is, for which we utilized equations (BIO) and (Bll) to derive 
(B12) and (B13) . Now, though, we require that subdivided droplets have a 
constant average second core volume moment fraction f 2 , 


— e /n + 2(u /n ) (u./n.) + £ i /n i W 2 
^ 2_ _ 1 l x^ i J J 1 j i 
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where and W ? +1 are the core volume second moments assigned to the sub- 

divided droplets. Equation (B15) , combined with (Bll), results in a fixed 
relative core volume variance for the droplets that are fractionated from 
intermediate-size coagulated particles. This assumption leads directly to the 
second core volume moment coagulation production rates: 
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We also have the net core volume second moment loss rate for droplets of size 
i due to coagulation. 
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The core volume moment coagulation terms ((B12), (B13) , and (B16)-(B18)) 
can be put in a semi-implicit numerical form similar to that of the droplet 
coagulation terms in the subsection. The Numerical Solution for Particle 
Coagulation; this can be done so that core volume is conserved exactly, and 
so that the core volume product u i u j the second moment equations (B16) 

and (B17) is not implicitly coupled to the first moment equation (this allows 
a simultaneous finite-difference solution for n, u, and e) . 


For the coagulation of condensation nuclei with aerosol droplets, the 
related coagulation terms for u and e are similar to those already given, 
except that for nuclei, u^ = n^V^ and = n^V-^ 2 . Moreover, the coagulation 
of a nucleus with a droplet is assumed to form only a larger droplet. We will 
not write out all of the lengthy numerical coagulation expressions here. 
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DROPLET WATER VAPOR GROWTH 


Droplet cores have a definite role in limiting the extent of water vapor 
growth that occurs when droplets move between physical environments. Equa- 
tion (8) in the main body of the report indicates that a droplet with a fixed 
core size will move along a trajectory in v - z space such that 


, _ dv , 

<p 3 — <p 

v dz z 


(B19) 


where the cj) 1 s are particle fluxes. In the aerosol continuity equation, the 
calculation of the total particle flux divergence in v - z space leads to 
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(B20) 


where 4> now represents just the vertical particle flux, which is defined by 
equation (24) in the text (but now has the units, particles cm“ 2 sec -1 ym -3 
yin -3 ). For a given droplet core size, the volume divergence term in equa- 
tion (B20) , which describes water vapor growth, can be written as 
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Integrating equation (B21) over all the core sizes, we obtain: 

h 2 o 

growth 

where $ (u) is the vertical core volume flux (ym 3 cm -2 sec -1 yirT* 3 ). 
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The corresponding results for u and e are easily derived, and are: 
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where it is the third core volume moment. 
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Obviously, we require a closure relation for tt to solve these 
equations. In our model, we only require a value for e in the regions of 
the atmosphere where droplets are evaporating; that is, high in the strato- 
sphere where water-vapor growth effects are small (because the product pW 
is nearly constant with height there) . To be consistent with the core evapo- 
ration routine that is outlined in the following section, we achieve closure 
by assuming a log-normal core size distribution (see eq. (B26) below) which 
leads to an approximation for tt. 


TT 




(B25) 


Equation (B25) may be substituted into equation (B24) . 


DROPLET EVAPORATION 


When droplets evaporate they release their cores. In equation (45) we 
determine the rate at which cores of all sizes evaporate from droplets of 
size i; let us call this rate q^ (cores cm -3 sec -1 ). We assume that all 
cores regardless of size are released at the same rate; hence, the size dis- 
tribution of any collection of incrementally evaporated cores is the same as 
that of the cores left in the droplets. To describe the sizes of the cores 
inside the droplets of size i, we have chosen a particular log-normal proba- 
bility distribution. 
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where is a normalization constant, W^ is the average size of the cores 

in the droplets, and a ^ is the core volume variance given by 


2 Zn 2 o . 
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The cores that are evaporating from droplets of size i must be allotted 
to condensation nuclei volume bins, m < i. Our geometric size bins have equal 
log-volume spacing. Accordingly, we can estimate the relative probability 
that an evaporated core is of size m by using 
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We want to conserve the number and volume of cores exactly during evaporation, 
so we proceed as follows. We divide the accessible core bins into large and 
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small sizes using the average core volume as a demarcation point; that 

is, we determine in- where 
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m io +1 


(B29) 


Then we renormalize the probability distributions for the large and small core 
sizes independently. 
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and calculate the average large and small core volumes and Wg_^, respec- 

tively, as, for example, 
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These volumes satisfy the inequality. 
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Next, we define the volume fractions. 
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which we use with the renormalized probabilities to compute the conden- 

sation nuclei production rates: 

,u> _ 
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Finally, we sum over all evaporating droplet sizes i to obtain the total 
nuclei production rates in each size category m. 
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To test the accuracy of our core evaporation algorithm, we have made 
several model calculations in which coagulation is omitted, and the gravita- 
tional settling velocity is made independent of particle size. Ihder these 
conditions the nucleation, growth, sedimentation, diffusion, rainout, and 
evaporation of droplets should not affect the size distribution of the total 
ensemble of cores at each altitude; that is, if we take the droplets at any 
height and evaporate them, the resulting collection of cores, when combined 
with the surrounding condensation nuclei, should have the same size distribu- 
tion as the condensation nuclei at the ground. Accordingly, the predicted 
size spectrum of the condensation nuclei in the evaporation regime above the 
aerosol layer should reproduce the nuclei size distribution at the ground. 

The results of our calculations for several different assumed condensa- 
tion nuclei distributions at the ground are given in figure 18. The model 
predictions are obviously in excellent agreement with the expected values. In 
the monodispersed case, less than one core in 10 13 has leaked outside the ori- 
ginal 0.1-pm size bin. Curve u d" of figure 18 shows the effects of coagula- 
tion and differential sedimentation on the core size distribution. On the 
average, droplet coagulation leads to bigger cores in the small size range, 
while gravitational settling effectively depletes the cores in the large-size 
range. 
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APPENDIX C 


A SOLUTION OF THE AEROSOL FINITE DIFFERENCE EQUATIONS 


In this appendix we outline a numerical scheme to solve the finite dif- 
ference continuity equations (63), which describe the evolution of aerosol 
particle concentrations. Some general comments about this scheme are in 
order. The equations for n^ 1 (or u^ 1 or e^ 1 ) as they are developed in 
the text, contain only the implicit concentrations of the smaller size acid 
droplets and of the same size condensation nuclei. The equations for the 
condensation nuclei contain nucleation, coagulation, and transport terms, as 
well as a production term due to droplet evaporation. If we treat this latter 

term as a source of condensation nuclei, a source that is fixed at the begin- 
ning of a time step, we can decouple the nuclei equations from the droplet and 

core volume moment equations. In this case, the solution strategy at each 

time step is to first solve the condensation nuclei continuity equation and 
then solve, in the order of increasing particle sizes, the aerosol droplet 
and core volume moment equations simultaneously, utilizing the nuclei concen- 
trations already calculated to evaluate implicitly nucleation and coagulation 
terms. The rates of production of nuclei by droplet evaporation can then be 
computed for the start of the next time-step cycle using the implicit core 
volume moments, u^ +1 and , to determine the evaporated core size spectrum. 

An alternative scheme is a sequential solution of the aerosol continuity 
equations, in which the variables are solved in the order n-n-u-e . The 
sequential approach has an advantage in that it allows more symmetry to be 
used in the implicit finite difference terms, particularly those terms 
describing core coagulation and water vapor growth, but it has a disadvantage 
in that it is more lengthy and costly. 

The finite-difference equations for aerosol droplet concentrations (or 
for condensation nuclei concentrations or core volume moments) may all be 
written in the general form: 
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(where D^ is not to be confused with a diffusion coefficient). Our solu- 
tion of these equations is based on the fact that for a given size category 
i the system of equations (Cl) can be represented by a tridiagonal matrix 
coupling together the aerosol number densities at different altitude levels k. 
We solve for the aerosol number densities by inverting the matrix (e.g., 
ref. 122). First, we define a relationship which couples the concentrations 
at adjacent altitude levels. 
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Next, using equations (Cl) and (C2) , we derive the recursion formulas 
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The lower boundary conditions are applied to determine and F.q which, 

together with equations (C3) , are used to calculate the E and F values at 
the higher altitude levels, k > 1. The upper boundary conditions are combined 
with the computed values of E and F for that level to determine n^T 1 (where 
k^=^K denotes the uppermost model level). Finally, working downward, the 
n^ are calculated in descending order using equation (C2) and the known 
E and F coefficients. 


The lower boundary condition in our model is a fixed concentration of 
condensation nuclei in each size bin and a removal flux of aerosol droplets. 
In the former case, we have 



n . 
xo 
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where n. is the fixed nuclei abundance, 
xo 

For aerosol droplets, on the other hand, we define an escaping flux. 
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The first component of <j) . is a sedimentation flux, and the second is a sur- 
face recombination flux (nominally, v B = 1 cm sec -1 in our calculations). The 
flux given by equation (C5) can be substituted into the aerosol continuity 
equation (63) to derive the appropriate values of E^ and F^ for the drop- 
lets . 


At the upper boundary of the model, we specify zero fluxes for all of the 
particles. 

In our calculations we assume that there are no particles smaller than 
our smallest size bin, and hence there is no growth or coagulation into the 
smallest size bin. We also prevent particles from growing or coalescing out 
of the largest size bin, although this mass flow would be negligible anyway. 
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( 


The finite-difference aerosol equations are connected to the gaseous 
photochemical equations by H 2 SO 4 condensation and evaporation; the mathemati- 
cal connection is made through equations (34) and (35). The numerical analogs 
of the sulfuric acid production and loss rates in equation (35) are formally 
related to equations (43) and (46). A model calculation normally proceeds as 
follows: For each time step, the evolving H 2 SO 4 vapor concentrations are com- 

puted using our photochemistry module which includes sulfuric acid condensa- 
tion and evaporation rates based on the aerosol droplet concentrations at the 
start of the step. The new H 2 SO 4 concentrations are then used to calculate 
the aerosol droplet formation and growth rates, and the continuity equations 
for the aerosol droplets, condensation nuclei, and core volume moments are 
solved over the same time interval. The numerical cycle is repeated until the 
desired simulated time has elapsed. 
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TABLE 1.- REACTIONS OF SULFUR COMPOUNDS 



Reaction 

Rate coefficients * 3 

Reference 

R1 

S + O 2 -*• SO + 0 

2.2X10 " 12 

62 

R2 

so + o 2 -> so 2 + 0 

3.0*10 - 13 e - 280o/T 

63 

R3 

SO + O 3 ■** SO 2 + O 2 

2.5X10 - 12 e - 105 °/ T 

63 

R4 

SO + N0 2 -*■ S0 2 + NO 

1.5 X 10 -11 

63 

R5 

S0 2 + 0 + M -»■ S0 3 + M 

3.4X10 -32 e _1130 / T (M = N 2 ) 

64 fc 

R 6 

S0 2 + H 0 2 S 0 3 + OH 

1.0X10 " 13 e - 141 °/ T 

65 ° 

R7 

S0 2 + OH + M HS0 3 + M 

8.2xl0 - 13 / (7.0 X 10 17 + [M]) 
(M = N 2 , 220 K) 

66 d 

R 8 

HS0 3 + OH -> S0 3 + H 2 0 

lx 10" 11 

e 

R9 

S0 3 + H 2 0 H 2 S0 4 

9.1X10 ' 13 

15 

RIO 

OCS + 0 SO + CO 

3.0X10 " 11 e - 227o/T 

-P 

67- 

Rll 

SO 2 washout 

3. 8 X 10 -6 (1 - 3 1 ~ z) z < 13 km 

9 



0 z > 13 km 


R12 

HSO 3 -> washout 

Rl 2 = Rl3 

h 

R13 

H 2 SO 4 -+ washout 

2.3X10 -5 - Z ^ z < 13 km 

i 



0 z > 13 km 


Jl 

S0 2 + hv + SO + 0 

1.9C-13)* 7 ’, 3.1 (- 10 ), 

68 



2 - 3 (— 8 ) , 4.1(-7), 2 . 7 (- 6 ) , 
8.5 (- 6 ) , 1.8(-5), 2 . 6 (-5) , 
3.3 (-5) 

69 

J 2 

OCS + hv -* S + CO 

1.5(-14) J ’, 1.9 (-11) , 

70 j 



1.4(-9) , 3.0(-8) , 2 . 8 (-7) , 
1.4(-6) , 5.0 (- 6 ) , l.l(-5), 
1 • 7 (-5 ) 

71 

J3 

H 2 SO 4 + hv SO 2 + products 

8.5(-16) J ', 1.6(-12) , 
1.3(-10), 2.4 (-9) , 1.5(-8) , 
4. 4 (- 8 ) , 8 . 6 (- 8 ) , 1. 3(-7) , 
1.6(-7) 

k 


a-k 


See footnotes on following page. 
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FOOTNOTES FOR TABLE 1 


a Rate coefficient units are sec” 1 for unimolecular , cm 3 sec” 1 for 
bimolecular , and cm 6 sec"* 1 for termolecular processes. 

£* Atkinson and Pitts (ref. 72) have obtained a similar rate constant, 
9.2*10 ~ 32 e "( 1000 ± 20 °)/T C m 6 sec" 1 , for M = N^O, and Westenberg and 
deHaas (ref. 73), a value of l.lxlO ” 31 e ”^ 1400 ± 50 )/ T C m 6 sec " 1 for M = He. 

^The rate constant measured at room temperature is 9*10 ” 16 cm 3 sec" 1 ; 
we have assumed a prefactor of lxlO ” 13 cm 3 sec " 1 and the corresponding acti- 
vation energy. Reaction R 6 is much less important than reaction R7 in the 
aerosol layer. 

^The rate constant value is based on a review of the available measure- 
ments, principally those by Castleman et al. (ref. 15). The observed pres- 
sure dependence has been simulated by a two-state reaction mechanism. 

e An estimated value; complemented by R9, reaction R 8 is proposed to 
replace a complex, but as yet undetermined, photoreaction sequence starting 
with HSO 3 and presumably leading to sulfuric acid vapor production. 

/k rezenski et al. (ref. 74) have obtained a similar rate coefficient of 
1.6X1CT 11 e” 2250 /’ 1 ’ cm 3 sec -1 . 

^The assumed SO 2 washout rate is based on an estimated low-altitude 
residence time of about 3 days (refs. 41, 47, 53). 

^Equal washout rates for HSO3 and H2SO4 are assumed. 

^The H 2 SO 4 washout rate is based on an assumed 0.5 day residence time 
near the surface. 

<7Twenty-four hour average photodissociation rates (sec" 1 ) are given for 
altitudes of 10, 15, 20, 25, 30, 35, 40, 45, and 50 km, respectively. The 
number 1.9(-13) = 1.9xl0~ 13 . 

^The H 2 SO 4 photoabsorption spectrum is assumed to be the same as the 
HC£ spectrum, with cross section data taken from reference 75. Ultraviolet 
decomposition of H 2 SO 4 in the upper stratosphere is assumed to lead to S 0 2 
production. 
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TABLE 2 .- A SUMMARY OF THE MAJOR AERONOMIC ASSUMPTIONS MADE IN THIS STUDY 


Assumption (parameters) 

_ . a 

Parameter uncertainty 

Physical effects 

Fixed atmospheric tem- 
perature and density 
profiles (air tempera- 
ture and density) 

These parameters vary 
with season and geo- 
graphical location but 
can be changed in the 
model accordingly. 

Large perturbations 
such as stratospheric 
warmings may also occur. 

Changes in temperature 
can affect the H2SO4 
supersaturation ratio 
and hence the height of 
the top of the aerosol 
layer. For small 
excursions in tempera- 
ture and density, the 
effect on chemical 
reaction rates is 
negligible. 

Parameterized one- 
dimensional transport 
by diffusion (diffu- 
sion coefficient) 

About a factor of 2-3 

The height of the 
equivalent tropopause 
level in the model, 
which directly affects 
the predicted height of 
the aerosol layer, is 
defined largely by the 
diffusion profile mini- 
mum. The diffusion 
coefficient also 
affects the transport 
rate of S 0 2 and OCS to, 
and the removal rate of 
aerosol from, the 
stratosphere. 

Reaction of OH with SO2 
as the rate limiting 
step in SO2 oxidation 
(overall SO2 oxidation 
rate) 

Unknown, but probably 
unimportant 

If the H2SO4 formation 
time were very long 
(~ years, compared to 
weeks for the OH reac- 
tion) the distribution 
of the sulfur gases and 
the height and shape of 
the aerosol layer would 
be changed signifi- 
cantly. 


a The parameter uncertainty factors quoted in this table are our estimates 
based on our evaluation of information currently available. 
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TABLE 2 .- CONTINUED 


Assumption (parameters) 

Parameter uncertainty a 

Physical effects 

Fast reaction of SO3 
with H2O to form H2SO4 
(rate constant for 
S 0 3 + H 2 0 ) 

The observed rate con- 
stant is fast, but for 
this type of reaction 
it could conceivably be 
orders of magnitude 
slower. 

Probably small, unless 
the rate were very slow, 
<10~ 17 cm 3 sec" 1 . 

An H2SO4 photolysis 
rate equal to that of 
HCL (H2SO4 photodisso- 
ciation coefficient) 

An order of magnitude 

Little effect, except 
possibly at the top of 
the aerosol layer where 
H2SO4 evaporates from 
droplets. 

Omission of hetero- 
geneous chemical inter- 
actions, including the 
condensation of HSO3 on 
particles (heteroge- 
neous reaction rates) 

For the overall rate of 
SO2 oxidation, unknown. 
For the overall chemi- 
cal rates of destruc- 
tion of species like 
°3 , N 2 0 5 , and C 10 N 0 2 , 
a factor of 2 in the 
extreme case. 

Heterogeneous chemistry 
could speed SO2 oxida- 
tion, which would have 
little consequence for 
our model. Surface 
effects on other gas- 
eous constituents could 
be significant, but the 
influence on the aero- 
sol would be second 
order. 

Fixed vertical profiles 
of OH, H 0 2 , H 2 0 , N 0 2 , 
and O3 (species con- 
centrations) 

, 

Since the species pro- 
files are calculated 
with a model including 
aerosol interactions , 
only the normal uncer- 
tainties associated 
with these species need 
to be considered. The 
most important uncer- 
tainties are in OH 
(± factor of 3 ) and H2O 
(± factor of 2 ). After 
a volcano, the species 
profiles could be 
affected by the large 
amounts of emitted sul- 
fur gases. 

OH affects the rate of 
SO2 oxidation, but for 
its probable range of 
values the effect would 
be small. H2O partly 
controls the H2SO4 
supersaturation in air 
and it could affect the 
top and bottom of the 
aerosol layer. In a 
volcanic plume, OH 
could be severely 
depleted by reaction 
with SO2 , and the aero- 
sol formation rate 
slowed significantly. 
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TABLE 2.- CONCLUDED 


Assumption (parameters) 

„ . a 

Parameter uncertaxnty 

Physical effects 

Ground-level abundances 
of OCS and S0 2 (concen- 
trations of S0 2 and 
OCS) 

The surface OCS concen- 
tration is probably 
uncertain by about 
10-20% . S0 2 at the 

ground is uncertain by 
at least a factor of 3, 
and more so at the 
tropopause (there is 
only one preliminary 
observation) . 

OCS is responsible for 
a large fraction of the 
aerosol sulfate mass 
above the tropopause. 
The S0 2 contribution to 
the aerosol layer by 
vertical transport is 
roughly proportional to 
its concentration near 
the tropopause. 

Omission of CS 2 as a 
source of stratospheric 
sulfur (deposition rate 
of sulfur in the 
stratosphere) 

Unknown, but it proba- 
bly results in less 
than a 20% uncertainty 
in the total rate of 
sulfur input to the 
stratosphere. 

An additional source of 
stratospheric sulfur 
would increase the 
amount of sulfate mass 
in the aerosol layer. 
However, CS 2 concentra- 
tions in the upper 
troposphere and strato- 
sphere are undetermined. 

Washout lifetimes of 
about 3 days for S0 2 » 
0.5 day for H 2 S 04 , and 
1 day for aerosol par- 
ticles near the ground 
(washout rate) . 

At least a factor of 
2-3 

For H 2 S 04 and aerosol, 
washout is an important 
sink. Without it, the 
aerosol layer would be 
much denser near 
the tropopause. S0 2 
washout also controls 
its abundance at the 
tropopause. 
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TABLE 3 .- A SUMMARY OF THE MAJOR AEROSOL PHYSICAL ASSUMPTIONS MADE IN 

THIS STUDY 


Assumption (parameters) 

Parameter uncertainty^ 

Physical effects 

Omission of homogeneous 
nucleation (aerosol 
nucleation rate) 

The homogeneous nuclea- 
tion rate may be uncer- 
tain by several orders 
of magnitude for strato- 
spheric conditions. 

Previous work has shown 
that homogeneous nucle- 
ation, or nucleation 
onto ions, is probably 
not significant under 
normal stratospheric 
conditions. 

A fixed heterogeneous 
heteromolecular nucle- 
ation rate of 10 -6 sec” 1 
(heterogeneous nucle- 
ation rate) 

A factor of 10 or more 

Theory suggests that 
heterogeneous hetero- 
raolecular nucleation is 
fast. Our assumed 
rate, ± factor of 10, 
would have little 
effect on our results. 

H2SO4-H2O solution 
vapor pressures from 
reference 108 (H2SO4 
and H 2 0 equilibrium 
vapor pressures) 

An order of magnitude 
or more at some alti- 
tudes (based on the 
disagreement between 
theory and measure- 
ments) 

In the aerosol layer, 
H2SO4 vapor supersat- 
urations are large 
enough that a change by 
a factor of 10 in the 
i H2SO4 vapor pressure 
has little consequence. 
The H2SO4 vapor pres- 
sure does affect the 
predicted height extent 
of the aerosol layer, 
and the importance of 
homogeneous nucleation. 

Use of vapor pressure 
for pure H2SO4-H2O 
solutions (H2SO4 and 
H2O equilibrium vapor 
pressures) 

A factor of 10 or more 
decrease in the H2SO4 
vapor pressure might 
occur when certain 
materials are dissolved 
in solution 

See discussion immedi- 
ately above 


a The parameter uncertainty factors quoted in this table are estimates 
based on our evaluation of information currently available. 
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TABLE 3 .- CONTINUED 


Assumption (parameters) 

Parameter uncertainty a 

Physical effects 

Liquid, not frozen, 
stratospheric particles 
(aerosol phase) 

Stratospheric aerosols 
are most probably solu- 
tion or slurry droplets 

Vapor condensation onto 
crystal surfaces would 
lead to modified parti- 
cle growth rates. The 
related vapor pressures 
over the solid parti- 
cles would, likewise, 
be different. 

The condensation rate 
of H2SO4 vapor on drop- 
lets given by refer- 
ence 106 (aerosol 
growth rate, H2SO4 
sticking coefficient) 

At low altitudes, the 
growth rate is probably 
uncertain by less than 
50 %. However, the 
H2SO4 sticking coeffi- 
cient is unknown, and 
it may lie anywhere 
between 0.03 and 1.0, 
which makes high- 
altitude growth rates 
significantly more 
uncertain. 

If the H2SO4 condensa- 
tion rate is altered, 
the H 2 S0 4 vapor concen- 
tration adjusts itself 
to maintain about the 
same net particle 
growth rate. Particle 
evaporation at high 
altitudes is affected 
in the same way. 

Droplet cores are inac- 
tive and do not grow by 
absorption and reaction 
of trace gases (droplet 
core volume, solution 
reaction rates) 

Unknown 

Core growth can affect 
aerosol size distribu- 
tions; its treatment 
requires detailed infor- 
mation about reactants 
and reactions which is 
not yet available. 

Aerosol droplets are in 
local thermal equilib- 
rium at the ambient air 
temperature (particle 
temperature) 

Probably negligible 

For stratospheric 
particles a change in 
thermal energy — due to 
the latent heat of con- 
densation or evapora- 
tion, radiation absorp- 
tion, or a change in 
ambient temperature — 
should be rapidly 
equilibrated by heat 
transfer with the sur- 
rounding air. 
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TABLE 3.- CONTINUED 


Assumption (parameters) 

. CL 

Parameter uncertainty 

Physical effects 

Coagulation kernels 
from reference 111 
(coagulation kernels) 

May be as large as a 
factor of 2 for some 
size ranges, but prob- 
ably less than 20% 
overall — based on com- 
parisons with experi- 
mental data. Atmo- 
spheric parameters such 
as temperature, den- 
sity, and viscosity 
affect the kernel 
values . 

Coagulation is an 
important sink for 
small particles near 
the tropopause level. 
Coagulation affects the 
total particle mixing 
ratio profile. Coagu- 
lation is an important 
growth mechanism in 
some size ranges. 

Sedimentation veloci- 
ties based on the 
Stokes-Cunningham for- 
mula (sedimentation 
velocity) 

About 20% or less; the 
determining equation 
has been studied exten- 
sively and is not very 
sensitive to variabil- 
ity. Additional uncer- 
tainties are introduced 
with air density and 
viscosity data. 

Sedimentation controls, 
to a large extent, the 
particle residence 
times in the strato- 
sphere, especially at 
higher altitudes. 

An inverse fourth power 
(r™*^) Aitken nuclei 
size distribution at 
the ground (Aitken 
nuclei concentration) 

Probably an order of 
magnitude, particularly 
in the small and large 
size ranges. The r -i+ 
distribution is one of 
many which can be fit 
to observational data. 

Small Aitken nuclei are 
a major source of new 
aerosol particles in 
the stratosphere. 

Large Aitken nuclei 
contribute to the large 
particle concentration 
in the lower strato- 
sphere. 

One-dimensional parti- 
cle transport by diffu- 
sion and sedimentation 
(net transport flux) 

About a factor of 2. 

In addition, horizontal 
transport and inhomoge- 
neities, as well as 
large scale 
stratospheric- 
tropospheric exchange, 
are not simulated by 
a one-dimensional 
model. 

In a model, vertical 
particle diffusion is 
most important in the 
lower stratosphere, 
where the diffusive 
fluxes are most uncer- 
tain. Sedimentation 
and diffusion oppose 
each other at higher 
altitudes and establish 
the particle mixing 
ratios above the aero- 
sol layer. 
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TABLE 3.- CONCLUDED 


Assumption (parameters) 

a 

Parameter uncertainty 

Physical effects 

An aerosol washout rate 
of about 1 day' 1 near 
the ground and about 
1 week -1 near the 
tropopause (washout 
rate) 

Factor of 2-3 


Washout is an important 
sink for stratospheric 
particles reaching the 
troposphere. However, 
a factor of 2 change in 
washout rate has little 
effect on the sulfate 
layer except near the 
tropopause. 

Omission of particle 
electric charge effects 
on aerosol growth and 
coagulation processes 
(growth and coagulation 
rates, particle charge) 

This omission probably 
leads to very small 
uncertainties in the 
physical growth rates. 
Stratospheric particle 
charges are not well 
determined . 

For growth, the conden- 
sation rate of polar 
molecules and their 
corresponding vapor 
pressures could be 
slightly altered. For 
coagulation, coulombic 
forces would probably 
not affect particle 
trajectories except for 
very small particles. 
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TABLE 4.- A SUMMARY OF THE MAJOR NUMERICAL APPROXIMATIONS MADE IN THIS STUDY 


Approximation 

Numerical uncertainty 

Computational effects 

Doubling of particle 

Some numerical disper- 

Volume doubling gives a 

volume between model 

sion may be introduced 

crude particle size 

size bins 

into particle coagula- 
tion and growth rates 

resolution. Coagulat- 
ing and growing parti- 
cles must often be 
subdivided into two 
adjacent size bins. 

A semi-implicit inte- 

The major uncertainty 

The semi- imp lie it 

gration scheme 

arises in droplet 
growth and sedimenta- 
tion processes. In 
these cases individual 
particles may eventu- 
ally be dispersed over 
several grid points. 

For a smooth, continu- 
ous distribution of 
particles, however, the 
net spreading should 
not introduce a signif- 
icant error. 

scheme leads to greater 
solution stability, but 
enhances numerical dif- 
fusion between grid 
points . 

Boundary conditions: 

With the numerical 

Boundary conditions 

no mass flow through 

boundaries at 0 and 

help to establish the 

the model boundaries 

58 km altitude and at 

net rates of mass flow 

except at the ground 

about 0.009 and 2.9 pm 

for the model. They 

where Aitken nuclei 

radial size, most of 

also partially deter- 

enter and aerosol drop- 

the stratospheric aero- 

mine the distributions 

lets leave, and where 

sol is sufficiently 

of condensation nuclei 

sulfur gases originate 

isolated from the 
boundaries to prevent 
large uncertainties 
from propagating into 
the results. Other 
types of boundary con- 
dition uncertainties 
(e.g., Aitken nuclei 
concentrations and gas 
abundances) are dis- 
cussed in tables 2 
and 3. 

and gaseous sulfur con- 
stituents in the upper 
troposphere and lower 
stratosphere. 
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TABLE 4.- CONCLUDED 


Approximation 

Numerical uncertainty 

Computational effects 

Treatment of particle 
growth and coagulation 
using discrete particle 
sizes 

No more than 10 to 20% 
in growth and coagula- 
tion rates 

Some numerical disper- 
sion is introduced into 
the growth and coagula- 
tion rates. Otherwise 
there appears to be 
little effect on the 
model predictions. 

A scheme to suppress 
numerical diffusion in 
the aerosol sedimenta- 
tion and growth rates 

Less than 10 to 20% in 
predicted particle 
height and size 
distributions based on 
modeling experiments 

Artificial numerical 
diffusion affects the 
settling and growth 
rates of particles, 
especially the larger 
particles. 

The technique used to 
solve stiff, coupled, 
nonlinear gaseous 
species continuity 
equations. 

For steady-state cal- 
culations, less than 
5%; for time-dependent 
calculations, less than 
10%. 

See published discus- 
sions of our chemistry 
model for more details. 




TABLE 5.- PHYSICAL PROPERTIES OF A 0.3-ym RADIUS AQUEOUS SULFURIC ACID 

SOLUTION DROPLET IN THE ATMOSPHERE 


Altitude, 

km 

Air 

tempera- 

ture, 

K 

Water vapor 
concentration , 
molecules cm” 3 

h 2 so 4 -h 2 o 

solution 
density, 
g cm” 3 

H 2 S0 4 

equilibrium 

vapor 

concentration, 

o 

molecules cm 

h 2 so 4 

weight 
fraction 
in solution, 
percent 

0 

288 

1.6 

(17) 

1.39 

2.48 

(5) 

48 

2 

275 

5.8 

(16) 

1.41 

4.78 

(4) 

50 

4 

262 

2.0 

(16) 

1.42 

4.96 

(3) 

51 

6 

249 

6.8 

(15) 

1.42 

2.32 

(2) 

51 

8 

236 

1.9 

(15) 

1.42 

1.04 

(1) 

51 

10 

223 

5.2 

(14) 

1.40 

1.46 

(-1) 

49 

12 

219 

1.0 

(14) 

1.52 

7.98 

(0) 

61 

14 

216 

1.7 

(13) 

1.62 

1.93 

(2) 

69 

16 

216 

1.3 

(13) 

1.63 

3.78 

(2) 

70 

18 

216 

1.0 

(13) 

1.65 

6.79 

(2) 

72 

20 

217 

7.8 

(12) 

1.67 

1.98 

(3) 

73 

22 

219 

5.9 

(12) 

1.69 

9.54 

(3) 

75 

24 

221 

4.4 

(12) 

1.71 

4.17 

(4) 

77 

26 

223 

3.3 

(12) 

1.73 

1.63 

(5) 

79 

28 

225 

2.5 

(12) 

1.75 

5.79 

(5) 

80 

30 

227 

1.9 

(12) 

1.76 

1.83 

(6) 

81 

32 

230 

1.4 

(12) 

1.78 

8.22 

(6) 

83 

34 

234 

1.1 

(12) 

1.79 

4.37 

(7) 

84 

36 

239 

8.3 

(ID 

1.81 

2.81 

(8) 

86 

38 

245 

6.3 

(ID 

1.83 

1.88 

(9) 

89 

40 

250 

4.8 

(ID 

1.84 

7.64 

(9) 

91 
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TABLE 6.- CONDENSATION (GROWTH) KERNELS, g , FOR H 2 S0 4 -H 2 0 SOLUTION 

DROPLETS IN AIR 


Particle 
radii, 
lti-\ym 
ude, 
km 


Condensation kernels 

(in units of 10 -1 5 ym cm 3 molecule -1 sec -1 ) 


0.01 


0.02 


0.04 


0.08 


0.16 


0.32 

0.64 

1.28 

2.56 

5.99 

3.49 

1.88 

0.976 

6.29 

3.77 

2.07 

1.08 

6.79 

4.22 

2.36 

1.25 

7.55 

4.90 

2.82 

1.51 

8.17 

5.56 

3.31 

1.81 

9.17 

6.58 

4.08 

2.28 

7.40 

5.64 

3.70 

2.14 

6.58 

5.31 

3.75 

2.29 

6.76 

5.80 

4.36 

2.82 

6.86 

6.15 

4.93 

3.41 

6.84 

6.33 

5.38 

3.98 

6.72 

6.37 

5.67 

4.50 

6.59 

6.35 

5.85 

4.93 

6.47 

6.31 

5.96 

5.26 

6.35 

6.24 

6.00 

5.49 

6.26 

6.18 

6.02 

5.65 

6.15 

6.09 

5.98 

5.72 

6.04 

6.01 

5.92 

5.75 

5.92 

5.90 

5.84 

5.72 

5.78 

5.76 

5.72 

5.64 

5 . 66 

5.65 

5.62 

5.56 


0 

2 

4 

6 

8 

10 

12 

14 

16 

18 

20 

22 

24 

26 

28 

30 

32 

34 

36 

38 

40 


15.0 
13.8 
13.2 

13.1 
12.7 
13.0 

9.59 

7.90 

7.68 

7.50 

7.28 

7.02 

6.80 

6.61 

6.45 

6.33 

6.20 

6.08 

5.95 

5.80 

5.67 


14.5 

13.5 
13.0 
12.9 

12.6 

12.9 

9.53 

7.87 

7.66 
7.48 
7.27 
7.01 

6.79 
6.61 
6.45 
6.33 
6.20 
6.08 
5.95 

5.80 

5.67 


13.6 
12.8 
12.5 
12.5 
12.3 

12.7 
9.41 

7.79 
7.61 
7.45 
7.24 
6.99 
6.78 
6.60 
6 . 64 
6.33 
6.19 
6.08 
5.95. 

5.80 
5.67 


11.9 

11.5 
11.4 
11.7 

11.6 
12.2 

9.13 

7.63 

7.50 

7.37 

7.19 
6.96 
6.75 
6.58 
6.43 
6.32 

6.19 
6.07 
5.94 
5.80 
5.67 


9.09 

9.15 

9.45 

10.0 

10.3 

11.1 

8.55 

7.29 
7.26 
7.21 
7.08 
6.88 
6.70 
6.54 
6.41 

6.30 
6.17 
6.06 
5.94 
5.79 
5.67 


Values are given for every third particle size category in the model. 


TABLE 7.- COAGULATION KERNELS FOR A SULFURIC ACID AEROSOL AT 16 km ALTITUDE 
(The air temperature is 216 K and the air number density is 
3.5xl0 18 molecules cm” 3 ) 


Droplet 

radii, 

a 

Coagulation kernel, cm 3 particle 1 sec 1 

ym 

0.01 

0.02 

0.04 

0.08 

0.16 

0.32 

0.64 

1.28 

2.56 

0.01 

1.87 (-9) & 

3.14 (-9) 

8.14 (-9) 

2.50 (-8) 

7.88 (-8) 

2.29 (-7) 

5.82 (-7) 

1.32 (-6) 

2.80 (-6) 

0.02 


2.60 (-9) 

4.24 (-9) 

1.02 (-8) 

2.72 (-8) 

6.92 (-8) 

1.60 (-7) 

3.47 (-7) 

EttEKsM 

0.04 



3.33 (-9) 

4.84 (-9) 

9.71 (-9) 

2.08 (-8) 

4.41 (-8) 

9.13 (-8) 

1.86 (-7) 

0.08 




3.33 (-9) 

4.03 (-9) 

6.73 (-9) 

1.27 (-8) 

2.48 (-8) 

4.92 (-8) 

0.16 





2.50 (-9) 

2.69 (-9) 

4.08 (-9) 

7.24 (-9) 

1.37 (-8) 

0.32 






1.64 (-9) 

1.70 (-9) 

2.48 (-9) 

4.27 (-9) 

0.64 







1.10 (-9) 

1.15 (-9) 

1.64 (-9) 

1.28 








8.18 (-10) 

8.77 (-10) 

2.56 









6.86 (-10) 


^Values are given for every third 
°1.87 (-9) = 1.87xl0 -9 . 

The kernels are symmetric; i.e., 


particle 


K(r x ,r 2 ) 


size category in the model. 
= K(r 2 ,ri). 






TABLE 8.- TERMINAL SEDIMENTATION VELOCITIES FOR ATMOSPHERIC AEROSOL PARTICLES 


\ Particle 
\ radii, 
\um 
Alti\ 
tude, 






Sedimentation velocity, cm 

sec 

a 

L 






km \ 

0 . 

01 

0 . 

02 

0.04 j 

0 . 

08 

0 . 

16 

0 . 

32 

0 . 

64 

1 . 

28 

2. 

56 

0 

2.39 

(-5) b 

4.41 

(-5) 

9.57 

(-5) 

2.35 

(-4) 

6.68 

(-4) 

2.13 

(-3) 

7.83 

(-3) 

12.95 

(-2) 

1.30 

(-1) 

2 

2.99 

(-5) 

5.52 

(-5) 

1.18 

(-4) 

2.80 

(-4) 

7.63 

(-4) 

2.41 

(-3) 

8.45 

(-3) 

3.15 

(-2) 

1.38 

(-1) 

4 

3.77 

(-5) 

6.95 

(-5) 

1.46 

(-4) 

3.37 

(-4) 

8.79 

(-4) 

2.67 

(-3) 

9.12 

(-3) 

3.36 

(-2) 

1.47 

(-1) 

6 

4.85 

(-5) 

8.87 

(-5) 

1.83 

(-4) 

4.11 

(-4) 

1.03 

(-3) 

2.97 

(-3) 

9.85 

(-3) 

3.57 

(-2) 

1.56 

(-D 

8 

6.33 

(-5) 

1.16 

(-4) 

2.36 

(-4) 

5.15 

(-4) 

1.24 

(-3) 

3.41 

(-3) 

1.08 

(-2) 

3.85 

(-2) 

1.66 

(-1) 

10 

8.37 

(-5) 

1.52 

(-4) 

3.05 

(-4) 

6.51 

(-4) 

1.50 

(-3) 

3.94 

(-3) 

1.20 

(-2) 

4.15 

(-2) 

1.78 

(-1) 

12 

1.12 

(-4) 

2.15 

(-4) 

4.34 

(-4) 

9.12 

(-4) 

2.04 

(-3) 

5.07 

(-3) 

1.45 

(-2) 

4.81 

(-2) 

1.87 

(-1) 

14 

1.55 

(-4) 

3.08 

(-4) 

6.25 

(-4) 

1.30 

(-3) 

2.82 

(-3) 

6.65 

(-3) 

1.78 

(-2) 

5.58 

(-2) 

1.98 

(-1) 

16 

2.11 

(-4) 

4.21 

(-4) 

8.51 

(-4) 

1.75 

(-3) 

3.73 

(-3) 

8.42 

(-3) 

2.12 

(-2) 

6.21 

(-2) 

2.09 

(-1) 

18 

2.88 

(-4) 

5.76 

(-4) 

1.16 

(-3) 

2.37 

(-3) 

4.96 

(-3) 

1.09 

(-2) 

2.59 

(-2) 

7.07 

(-2) 

2.24 

(-1) 

20 

3.92 

(-4) 

7.86 

(-4) 

1.59 

(-3) 

3.22 

(-3) 

6 . 66 

(-3) 

1.48 

(-2) 

3.25 

(-2) 

8.29 

(-2) 

2.45 

(-1) 

22 

5.37 

(-4) 

1.08 

(-3) 

2.19 

(-3) 

4.43 

(-3) 

9.08 

(-3) 

1.91 

(-2) 

4.21 

(-2) 

1.01 

(-1) 

2.75 

(-1) 

24 

7.30 

(-4) 

1.49 

(-3) 

3.00 

(-3) 

6.07 

(-3) 

1.24 

(-2) 

2.56 

(-2) 

5.51 

(-2) 

1.26 

(-1) 

3.19 

(-1) 

26 

9.93 

(-4) 

2.04 

(-3) 

4.11 

(-3) 

8.30 

(-3) 

1.68 

(-2) 

3.46 

(-2) 

7.29 

(-2) 

1.61 

(-1) 

3.84 

(-1) 

28 

1.35 

(-3) 

2.78 

(-3) 

5.62 

(-3) 

1.13 

(-2) 

2.29 

(-2) 

4.67 

(-2) 

9.71 

(-2) 

2.09 

(-1) 

4.78 

(-1) 

30 

1.82 

(-3) 

3.78 

(-3) 

7.63 

(-3) 

1.54 

(-2) 

3.10 

(-2) 

6.30 

(-2) 

1.30 

(-1) 

2.74 

(-1) 

6.07 

(-1) 

32 

2.44 

(-3) 

5.07 

(-3) 

1.03 

(-2) 

2.08 

(-2) 

4.19 

(-2) 

8.48 

(-2) 

1.73 

(-1) 

3.61 

(-1) 

7.79 

(-1) 

34 

3.29 

(-3) 

6.59 

(-3) 

1.37 

(-2) 

2.82 

(-2) 

5.69 

(-2) 

1.15 

(-1) 

2.33 

(-1) 

4.80 

(-1) 

1.01 


36 

4.41 

(-3) 

8.83 

(-3) 

1.77 

(-2) 

3.78 

(-2) 

7.67 

(-2) 

1.55 

(-1) 

3.13 

(-1) 

6.39 

(-1) 

1.33 


38 

5.83 

(-3) 

1.17 

(-2) 

2.33 

(-2) 

4.95 

(-2) 

1.02 

(-1) 

2.05 

(-1) 

4.15 

(-1) 

8.43 

(-1) 

1.73 


40 

7.72 

(-3) 

1.54 

(-2) 

3.09 

(-2) 

6.43 

(-2) 

1.35 

(-1) 

2.73 

(-1) 

5.51 

(-1) 

1.11 

(-1) 

2.27 



^Values are given for every third particle size category in the model 
z . 39 (-5) = 2.39xl0 -5 . 






Figure 1.- Carbonyl sulfide photoabsorption cross sections. Below 226 nm 
the data are taken from reference 70 at 232 K, and above 226 nm 
from reference 71 at 300 K. 



Figure 2.- Gas phase sulfur interactions that are included in our model. 
The dashed lines indicate possible heterogeneous reactions which 
are not treated. 
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Figure 3.- Eddy diffusion coefficients in the lower atmosphere. 


io 5 io 6 io 7 io 8 io 9 io 10 



Figure 4.- Neutral species concentrations up to 58 km. 
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MIXING RATIO, ppbv 


Figure 5.- Predicted sulfur compound mixing ratios in the atmosphere. 



Figure 6.- A partial schematic outline of the- physical and chemical 

interactions included in our one-dimensional model. C = coagulation, 
Ch = chemistry, D = diffusion, E = evaporation, Em = emission, 

G = condensation and growth, I = injection, N = nucleation, 

P = photolysis, S = sedimentation, W = washout and rainout. 
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Figure 7 
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Figure 8.- Aerosol droplet growth. The evolution of a droplet size 
distribution following a change in the droplet growth rate is shown 
Both analytical and numerical solutions are given. The conditions 
for the calculation are described in the text. 



Figure 9.- Aerosol coagulation. Calculated transient decay curves are 
shown for the number of particles in the smallest model size bin, 
and for the total number of particles of all sizes, when one particle 
is initially placed in the smallest model bin. Time is measured 
relative to the characteristic coagulation time, T . The exact 
solutions of Smoluchowski are shown for comparison. 

PARTICLE RADIUS, fJi m 


0.01 OJ i.O 3.0 



Figure 10.- The development of an aerosol size spectrum by coagulation. 
The conditions for the calculation are the same as for figure 9. 
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Figure 11.- Particle sedimentation and diffusion. The steady-state altitude 
distributions of 0.1 and 1.0 ym radius particles resulting from their 
continuous injection at 20 km are presented. Particle mixing ratios 
calculated for two different diffusion coefficients are shown. Theoreti- 
cal and numerical solutions are both given for comparison. See the text 
for details. 



Figure 12.- Particle diffusion. The time-dependent diffusive transport of 
aerosol particles released at 20 km is shown. The initial particle 
concentration is 100 cm" 3 . The diffusion coefficient is given in 
figure 3. 
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Figure 13.- Observed and calculated aerosol sulfate mixing ratios (yg SO^/kg 
air). The observations (ref. 14) were made with filters, and have been 
averaged over latitude from pole to pole. The error bars give the 
extreme values found at individual latitudes. 



10 100 1000 10,000 
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Figure 14.- Observed and calculated total particle number mixing ratios 
(particles/mg air). The data were obtained over a period of 15 years 
by various investigators using condensation chambers. The large 
spread in the data may be partly due to real differences in aerosol 
abundances and partly due to different sampling techniques. 
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LARGE PARTICLE (r > 0.15 jum) MIXING RATIO, PARTICLES/mg AIR 


Figure 15.- Observed and calculated mixing ratios of particles with radii 
>0.15 ym. The observations (ref. 9) were made using a light-scattering 
particle counter. A large number of measurements were made over 
Wyoming; the mean value and a measure of the standard deviation are 
presented for different seasons. 
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RATIO OF PARTICLES > 0.15 fin 1 RADIUS 
TO THOSE > 0.25 /xm RADIUS 


Figure 16.- Observed and measured concentration ratios of particles with 

radii >0.15 ym to those with radii >0.25 pm. The observations (ref. 119) 
were made with a light-scattering particle counter. Numerous flights 
were made over Wyoming and worldwide. The average value and a measure 
of the standard deviation are illustrated. 
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Figure 17.- Observed and calculated aerosol size distributions. The 
calculated distributions are for altitudes of 16 and 20 km. The 
observations are represented by three analytic functions that are 
based upon close fits to a variety of measured distributions 
(refs. 102 and 119). An r -i+ size distribution, typical of tropospheric 
aerosol, is also shown for comparison. 
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Figure 18.- Calculated condensation nuclei size spectra at 40 km altitude 
(above the model droplet evaporation level) . The assumed condensation 
nuclei distributions at the ground are: (a) monodispersed 0.1 pm radius 

particles, (b) an r” 3 size distribution, (c) and (d) an r“ 4 size distri- 
bution. The points within the symbols correspond to the numerical pre- 
dictions; the curves have been added for clarity. Results (a) to (c) 
were obtained using the aerosol model without droplet coagulation or 
differential sedimentation. Result (d) is for the complete aerosol model. 
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